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The dibaryon model for NN interaction, which implies the formation of an inter¬ 
mediate six-quark bag dressed by a cr-field, is applied to the 2>N system, where it 
results in a new three-body force of scalar nature between the six-quark bag and a 
third nucleon. A new multicomponent formalism is developed to describe three-body 
systems with nonstatic pairwise interactions and non-nucleonic degrees of freedom. 
Precise variational calculations of 2>N bound states are carried out in the dressed- 
bag model including the new scalar three-body force. The unified coupling constants 
and form factors for 2N and 2>N force operators are used in the present approach, 
in a sharp contrast to conventional meson-exchange models. It is shown that this 
three-body force gives at least half the 2>N total binding energy, while the weight of 
non-nucleonic components in the and ^He wavefunctions can exceed 10%. The 
new force model provides a very good description of 3A bound states with a reason¬ 
able magnitude of the aNN coupling constant. A new Coulomb 3A force between 
the third nucleon and dibaryon is found to be very important for a correct description 
of the Coulomb energy and r.m.s. charge radius in ^He. In view of the new results 
for Coulomb displacement energy obtained here for A = 3 nuclei, an explanation for 
the long-term Nolen-Schiffer paradox in nuclear physics is suggested. The role of the 
charge-symmetry-breaking effects in the nuclear force is discussed. 

1. INTRODUCTION. CURRENT PROBLEMS IN A CONSISTENT DESCRIPTION 
OE NN AND 3N SYSTEMS WITH TRADITIONAL EORCE MODELS 

A few historical remarks should be done at first. Current rather high activity in 
few-body physics started since the beginning of 1960-s, after mathematical formulation 
of the Eaddeev eqnations for three-body problem. The aim was claimed to establish 
nnambignonsly off-shell properties of the two-body f-matrix, which cannot be derived 
from two-body scattering data only. It has been hoped in that time that jnst accnrate 
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solving 3N scattering problem is able to put strong constraints for the off-shell properties 
of the two-nucleon t-matrix. However, more than forty years passed since that, but still 
now we are unable to formulate such a two-nucleon t-matrix, which can explain fully 
quantitatively the properties of even 3N systems. 

Moreover, from that time, many puzzles in few-nucleon scattering experiments have 
been revealed which could not be explained by the current force models based on Yukawa 
concept. Among all such puzzles, we mention here only the most remarkable ones, such 
as the Ay puzzle in iV -|- d and N + ^He scattering [1, 2], disagreements at the minima of 
differential cross sections (Sagara puzzle) at i? ~ 150-200 MeV and polarization data for 
N+d [3], N+d, N+d [4], and scattering, and many others. The strongest discrep¬ 

ancy between current theories and respective experiments has been found in studies of the 
short-range NN correlations in the ^He(e, e'pp) [5], "‘He( 7 ,pp) [6], and ^He(e, e'iViV) [7] 
processes. In addition to these particular problems, there are more fundamental problems 
in the current theory of nuclear forces, e.g., strong discrepancies between the nNN, nNA, 
and pNN form factors used both in one-boson-exchange (OBE) models for the description 
of elastic and inelastic scattering and in the consistent parametrization of 2N and 3N 
forces [8-11]. Many of these difficulties are attributed to a rather poor knowledge of the 
short-range behavior of nuclear forces. This behavior was traditionally associated with 
the vector cu-meson exchange. However, the characteristic range of this cu-exchange (for 
rriuj — 780 MeV) is equal to about — 0.2-0.3 fm, i.e., is deeply inside the internucleon 
overlap region. 

In fact, since Yukawa the nucleon-nucleon interaction is explained by a t-channel 
exchange of mesons between two nucleons. The very successful Bonn, Nijmegen, Argonne, 
and other modern NN potentials prove the success of this approach. But the short- and 
intermediate-range region in these potentials is more parametrized than parameter-free 
microscopically described. 

Besides of the evident difficulties with the description for short-range nuclear force 
there are also quite serious problems with consistent description of basic intermediate- 
range attraction between nucleons. In the traditional OBE models this attraction is 
described as a t-channel a-exchange with the artificially enhanced aNN vertices. How¬ 
ever, the accurate modern calculations of the intermediate-range NN interaction [12, 13] 
within the 27r-exchange model with the tttt s-wave interaction have revealed that this 
t-channel mechanism cannot give a strong intermediate-range attraction in the NN sec¬ 
tor, which is necessary for binding of a deuteron and fitting the NN phase shifts. This 
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conclusion has also been corroborated by recent independent calculations [14], Thus, the 
t-channel mechanism of the a-meson exchange should be replaced by some other alterna¬ 
tive mechanism, which should result in the strong intermediate-range attraction required 
by even existence of nuclei. 

When analyzing the deep reasons for all these failures, we must look on a most general 
element, which is common for all the numerous NN force models tested in few-nucleon 
calculations for last 40 years. This common element is just the Yukawa concept for 
the strong interaction of nucleons in nuclei. Hence, if, after more than 40 years of de¬ 
velopment, we are still unable to explain quantitatively and consistently even the basic 
properties of 3N and 4iV systems at low energies and relatively simple processes like 
pp — pp-y, this concept, which is a cornerstone of all building of nuclear physics, should 
be analyzed critically, especially in the regions where applicability of this concept looks 
rather questionable. 

Since the quark picture and QCD have been developed, the "nucleon-nucleon force 
community" is more and more convinced that at short ranges the quark degrees of freedom 
must play an important role. One of possible mechanisms for short-range NN interaction 
is the formation of the six-quark bag (dibaryon) in the s-channel. Qualitatively many 
would agree with this statement. But to obtain a quantitative description of the nucleon- 
nucleon and the few-nucleon experimental data with this approach with the same quality 
as the commonly used Bonn, Nijmegen, Argonne, and other equivalent potentials is a 
quite different problem. 

Within the 6q dynamics it has long been known [15-19] that the mixing of the com¬ 
pletely symmetric s®[6] component with the mixed-symmetry s^p^[42] component can 
determine the structure of the whole short-range interaction (in the A-wave)^^. Assum¬ 
ing a reasonable qq interaction model, many authors (see, e.g., [20-23]) have suggested 
that this mixture will result in both strong short-range repulsion (associated mainly with 
the component) and intermediate-range attraction (associated mainly with the above 
mixed-symmetry component). However, recent studies [22, 23] for NN scattering on 
the basis of the newly developed Goldstone-boson-exchange (GBE) qq interaction have 
resulted in a purely repulsive NN contributions from both s®[6] and s^p^[42] 6g compo¬ 
nents. There is no need to say that any quark-motivated model for the NN force with 

We will denote the NN partial waves by capital letters {S, P ..while the partial waves in all other 

cases will be denoted by small letters. 
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TT-exchange between quarks inevitably leads to the well-established Yukawa 7r-exchange 
interaction between nucleons at long distances. 

Trying to solve the above problems (and to understand more deeply the mechanism 
for the intermediate- and short-range NN interaction), the Moscow-Tiibingen group 
suggested to replace the conventional Yukawa meson-exchange (t-channel) mechanism (at 
intermediate and short ranges) by the contribution of a s-channel mechanism describing 
the formation of a dressed 6q bag in the intermediate state such as |s® -|- a) or -|- 
27r) [8, 24], It has been shown that, due to the change in the symmetry of the 6q state 
in the transition from the NN channel to the intermediate dressed-bag state, the strong 
scalar a-field arises around the symmetric 6g bag. This intensive a-field squeezes the bag 
and increases its density [25]. The high quark density in the symmetric 6g bag enhances 
the meson field fluctuations around the bag and thereby partially restores the chiral 
symmetry [26]. Therefore, the masses of constituent quarks and cr mesons decrease [8]. 
As a result of this phase transition, the dressed bag mass decreases considerably (i.e., a 
large gain in energy arises), which manifests itself as a strong effective attraction in the 
NN channel at intermediate distances. The contribution of the s-channel mechanism 
would generally be much larger due to resonance-like enhancement^). 

In our previous works [8, 24] on the basis of the above arguments we proposed a new 
dibaryon model for the NN interaction (referred further to as the “dressed-bag model” 
(DBM)),which provided a quite good description of both NN phase shifts up to 1 GeV 
and the deuteron structure. The developed model includes the conventional t-channel 
contributions (Yukawa tt- and 27r-exchanges) at long and intermediate distances and the 
s-channel contributions due to the formation of intermediate dressed-bag states at short 
distances. The most important distinction of such an approach from conventional models 
for nuclear forces is the explicit appearance of a non-nucleonic component in the total 
wavefunction of the system, which necessarily implies the presence of new three-body 
forces (3BF) of several kinds in the 3N system. These new 3BF differ from convention¬ 
ally used models for three-body forces. One important aspect of the novel 3BF should be 
emphasized here. In conventional OBE models, the main contribution to NN attraction 
is due to the t-channel a-exchange. However, the 3BF models suggested until now (such 

In the theory of nuclear reactions, the t-channel mechanism can be associated with the direct nuclear 
reaction, where only a few degrees of freedom are important, while the s-channel mechanism can be 
associated with resonance-like (or compound-nucleus-like) nuclear reactions with much larger cross 
sections at low energies. 
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as Urbana-Illinois or Tucson-Melbourne) are mainly based on the 27r-exchange with in¬ 
termediate A-isobar production, and the a-exchange is either not taken into account at 
all, or is of little importance in these models. In contrast, the a-exchange in our approach 
dominates in both NN and 3N forces. In fact, in our approach just the unified strong 
(T-field couples both two- and three (and more)-nucleon systems, i.e., the general pattern 
of the nuclear interaction appears to be more consistent. 

Our recent considerations have revealed that this dibaryon mode is extremely useful 
in the explanation of very numerous facts and experimental results in nuclear physics, in 
general. We note here only a few of them. 

1. The presence of dibaryon degree of freedom (DDF) can result in very natural ex¬ 
planation of cumulative effects (e.g., the production of cumulative particles in high- 
energy collisions [27]). 

2. DDF leads to automatic enhancement of near-threshold cross sections for one- and 
two-meson production in pp, pd, etc. collisions, which is required by many mod¬ 
ern experiments (e.g., the so-called ABC puzzle [28]). This is due to an effective 
enhancement of meson-dibaryon coupling as compared to meson-nucleon coupling. 

3. The incorporation of DDF makes it possible (without the artificial enhancement 
of meson-nucleon form factors) to share the large momentum of an incident probe 
(e.g., high-energy photon) among other nucleons in the target nucleus. 

4. The DDF produces in a very natural way a new short-range currents required by 
almost all experiments associated with high momentum and energy transfers. 

5. Presence of the dressed 6g bag components in nuclear wave functions leads automat¬ 
ically to a smooth matching between the nucleonic (at low momentum transferred) 
and quark currents (at very high momentum transferred) and, at the same time, 
results in a correct counting rules at high momentum transferred. 

So, it should be very important to test the above dibaryon concept of nuclear force 
in a concise and consistent 3iV calculations and to compare the predictions of the new 
model with the results of the conventional meson-exchange models. 

Thus, the aim of this work is to make a comprehensive study of the properties of 
the 3iV system with NN and 3N forces given by the DBM. However, DBM introduces 
explicitly non-nucleonic (quark-meson) channels. Therefore, it is necessary to introduce 
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a selfconsistent multichannel few-body formalism for the study of 3N system with DBM 
interaction. We develope in this work such a general formalism, based on the approach 
which was suggested in 1980-s by Merkuriev’s group [29-31] for the boundary-condition- 
type model for pairwise interactions. This general formalism leads immediately to a 
replacement of all two-body forces related to the dibaryon mechanism to the respective 
three (and many)-body forces, leaving two-body character only for long-range Yukawa 
vr and 27r exchanges, which are of little importance for the nuclear binding. Another 
straightforward sequence of the formalism developed here is a strong energy dependence 
of these many-body forces. In the work we study all these aspects in detail when applying 
to the 3N system properties. The preliminary version of this work is published in [32]. 

This paper is organized as follows. In Section 2, we present a new general multichannel 
formalism for description of two- and three-body system of particles having inner degrees 
of freedom. In Section 3, we give a brief description of the DBM for the NN system. 
In Section 4, we treat the 3N system with DBM interactions, including a new 3BF. In 
Section 5, some details of our variational method are discussed, including calculation of 
the matrix elements for new Coulomb 3BF. The results of our calculations for ground 
states of and ^He are given in Section 6 while in Section 7 we discuss the role of the 
new three-body force and present a new explanation for the Coulomb displacement energy 
in ^He within our interaction model. A comprehensive discussion of the most important 
results found in the work is given in Section 8. In the Conclusion we summarize the main 
results of the work. In the Appendix we give the formulas for the matrix elements of all 
DBM interactions taken in the Gaussian symmetrized variational basis. 

2. THE GENERAL MULTIGOMPONENT EORMALISM EOR 2N AND 3N 
SYSTEMS WITH GOUPLED INTERNAL AND EXTERNAL GHANNELS 

In 1980-s, the former Leningrad group (Merkuriev, Motovilov, Makarov, Pavlov, et al.) 
has constructed and substantiated with mathematical rigor a model of strong interaction 
with the coupling of external and internal channels [29-31, 33]. This model was a partic¬ 
ular realization of a general approach to interaction of particles having inner degrees of 
freedom. The basic physical hypothesis is that an energy-dependent interaction appears 
as a result of internal structure of interacting particles"^). A general scheme proposed by 


"‘I From more general point of view, the explicit energy dependence of interaction reflects its nonlocality 
in the time, while this time nonlocality, in its turn, is a result of some excluded degrees of freedom. So, 
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Merkuriev et al. has been based on assumption on existence of two independent channels: 
external one, which describes the motion of particles considered as elementary bodies, i.e., 
neglecting their inner structure, and internal one, which describes the dynamics of inner 
degrees of freedom. These channels can have quite different physical and mathematical 
nature and their dynamics are governed by independent Hamiltonians. The main issues 
here were - how to define the coupling between external and internal channels and how to 
derive corresponding dynamical equations (of Schrodinger or Faddeev type) for particle 
motion in the external channel. 

In [29-31] this coupling has been postulated via boundary conditions on some hyper¬ 
surface. Thus, such an approach is well applicable to hybrid models for NN interaction, 
which were rather popular in 1980-s, e.g., the quark compound bag (QCB) model sug¬ 
gested by Simonov [34]. As for the 3N system, the formalism of incorporation of the 
internal channels (6g bags) has been proposed for the first time also within the QCB 
model [35]. The general scheme by Merkuriev et al. has allowed to substantiate this 
formalism. 

In QCB-like models the coupling between the external (NN) and the inner (bag) chan¬ 
nels is given just on some hypersurface, similarly to the well-known i?-matrix approach in 
nuclear physics. Later on, such a general approach has been applied to the two-channel 
Hamiltonian model, where the internal Hamiltonian had pure discrete spectrum and the 
only restriction imposed on the operators coupling the external and internal channels 
was their boundness [33]. The above general multichannel scheme has straightforwardly 
been extended to three-body problem. In particular, it has been shown for the two above 
models that elimination of the internal channels leads to the following recipe for em¬ 
bedding the energy-dependent pair interactions into three-body problem: replacement of 
pair energy by difference between three-body energy and kinetic energy of third particle: 
Ea ^ E — ta [31, 33]. It has also been proved that the resulted Faddeev equations for 
external channel belong with such energy-dependent potentials to the Fredholm class and 
are equivalent to four-channel Schrodinger equation. 

Our aim here is to extend our new NN force model - DBM - by using the above 
Merkuriev et al. approach to the 3N system. There are external (nucleon-nucleon) 
and internal (quark-meson) channels in our model, and coupling between them is deter¬ 
mined within a microscopical quark-meson approach. In this section we present a general 

the explicit energy dependence is signalling about some inner hidden (e.g., quark) degrees of freedom 

in NN interaction. 
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multicomponent formalism for description of systems of two and three particles having 
internal structure, without assuming any specific form for coupling between the external 
and internal channels. 


2.1. Two-body system 


We adopt that the total dynamics in two-body system is governed by a selfconjugated 
Hamiltonian h acting in the orthogonal sum of spaces: 

where is the external Hilbert space of states describing motion of particles neglecting 
their internal structure and is the internal Hilbert space corresponding to internal 
degrees of freedom. Thus the total state of the system T G Ti can be written as a two- 
component column: 

/ ^ex 

T= . 

I 

The two spaces, Ti®"" and 7-f“, can have quite different nature, e.g., in the case of NN 
system depends on the relative coordinate (or momentum) of two nucleons and their 
spins, while can depend on quark and meson variables. The two independent Hamil¬ 
tonians are defined in each of these spaces: h®'' acts in Ti®"" and h™ acts in TY™. Here h®"" 
includes the kinetic energy of relative motion and some part of two-body interaction n®^: 

h^^ = t + n"". 



For NN system includes the peripheral part of meson-exchange potential and Coulomb 
interaction between nucleons (if they are protons). Coupling between external and in¬ 
ternal channels is determined formally by some transition operators: = (h™-®^)*. 

Further, one can write down the total Hamiltonian h as a matrix operator: 


h = 


/^ex 

h'n.ex 



( 1 ) 


not specifying so far the coupling operators (if operators and h™ are self-adjoint and 
/^ex,in jg bounded then the Hamiltonian h is the self-adjoint operator in Ti). 

Thus one can write down the two-component Schrodinger equation 


hT = ET, 
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and by excluding the internal channel wavefunction one obtains an effective Schrodinger 
equation in the external channel 

^ (^2) 

with an effective “pseudo-Hamiltonian”: 

= h"" + = t + v^^ + w{E), (3) 

which depends on energy E via the resolvent of internal Hamiltonian = (-E — 

(From mathematical point of view, an operator depending on the spectral parameter is 
not operator in all, because its domain depends on the spectral parameter. Thus, this 
object should not be called Hamiltonian. However, physicists do not turn their attention 
to the fact and use energy-dependent interactions very widely.) 

Having the solution of effective equation (2), one can “restore” the excluded internal 
state unambiguously: 

(4) 

2.2. Three-body system 

In three-body system we have three different internal spaces = 1,2,3) and one 

common external space The three-body internal space Ti™ is a direct product of 

the two-body internal space related to the pair (jk) and single-particle space describing 
motion of third particle (i). Here we use the conventional numbering of particles: {ijk) = 
(123), (231), (312). The own three-body Hamiltonian acts in each internal space as: 

Hr = h% ® I, + ® U; (ijk) = (123), (231), (312), (5) 

where EJl is the two-body internal Hamiltonian for the pair (jk), I is unity operator and 
ti is kinetic energy of third particle (i) in respect to the center of mass of the pair (jk). 
(Here and below we use capital letters for three-body quantities and small letters for 
two-body ones.) 

The external three-body Hamiltonian acts in the external space TCf and includes the 
total kinetic energy T and the sum of external two-body interactions, which were incor¬ 
porated to the external two-body Hamiltonians: 

nr = r + 5^ «•“. 

i<j 
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A state in the full three-body Hilbert space 

i 

can be written as a four-component column: 

/ \ 


'I /3 = 




) 

Thus, the total Hamiltonian, i^s, of the three-body system acting in H 3 can be written 
as (4 X 4) matrix: 

f /fex jjex,in jjex,in jjex,m ^ 

Hm,ex ^in q q 


H. = 


H^n,ex Q 
0 


0 


0 


( 6 ) 


Ht / 


Here we suppose: 

(i) there is no direct coupling between different internal channels TY™ and TY™ for i 7 ^ j; 

(ii) the channel coupling operators do not involve the third (free) particle: 

= (7) 


Writing the four-component Schrodinger equation with Hamiltonian ( 6 ): 

= -ET3, (8) 

and excluding three internal channels from it (it is simple due to the supposed absence of 
direct coupling between different internal channels), one obtains an effective Schrodinger 
equation for the external three-body wavefunction 

(9) 

with an effective (pseudo)Hamiltonian: 

Hf(E) = //r + E GfiE) (10) 

i 

where the resolvent of internal Hamiltonian G™ is a convolution of the two-body internal 
resolvent g^k of pair {jk) and the free motion resolvent for the third particle (i): 

00 

Gf = (E- //■”)-■ = T I = j>»(i5 _ (.) ® E.. (11) 
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Thus the effective Hamiltonian in external three-body channel takes the form: 

//" = r + + w,t(E - «,)}. (12) 

jk 

i.e., the total effective interaction in external channel of the three-body system is a sum of 
the two-body external potentials and the two-body effective interactions with replace¬ 
ment of pair energy Ei with difference between the total three-body energy and operator 
for the relative-motion kinetic energy of third particle: ^ E — ti. 

Just this recipe for inclusion of pair energy-dependent interactions in three-body prob¬ 
lem is widely used in Faddeev calculations. This recipe has been rigorously proved in 
works of Merkuriev et al. for two-channel model without continuous spectrum in internal 
channel [33] and, in particular, for the boundary condition model [31]. We see, however, 
that this result is a direct consequence of the above two assumptions and by no means is 
related to usage of any specific interaction model®^. 

The resulted form of the effective three-body Hamiltonian (12) is suitable for the 
Faddeev reduction. However, it should be emphasized that each term in the effective 
Hamiltonian (12) includes a dependence on the kinetic energy of the third particle, i.e., the 
each term is, generally speaking, a three-body force. In spite of three-body character 
of such effective potentials, the corresponding Faddeev equations have the Fredholm 
property and are equivalent to four-channel Schrodinger equation (it has been proved for 
model with discrete internal spectrum [33]). 

2.3. A new three-body force in the three-body system with external and internal 

channels 

In each internal channel one can introduce a new interaction between third particle 
and the pair as a whole. This leads to replacement of the operator for kinetic energy of 
the third particle ti by some (single-particle) Hamiltonian hy. 


ti ^ hi = ti + Vi (13) 


In the literature, however, there were also discussions of the alternative variants for embedding energy- 
dependent pairwise force into three-body system [36, 37]. These schemes suppose that the effective 
total energy £12 of the two-body subsystem in the three-body system is obtained from the total three- 
body energy E in the following way: £12 = E — ti 2 — (uis) — (^ 23 ), where Vij is the two-body interaction 
between particles i and j and an averaging is supposed with the exact 3N wavefunction. 
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in Eq. (5) for Hf^, viz.: 

= + (14) 

Physical meaning of snch interactions will be discnssed below and here we treat only the 
formal aspects of their introduction. As the internal Hamiltonian (14) is still a direct 
sum of the two-body internal Hamiltonian and the Hamiltonian corresponding to relative 
motion of the third nucleon, then its resolvent can be expressed as a convolution of two 
subresolvents: 

OO 

Gf = (B - Hr)-' = ^ / a%{^)3i(E - z)dz, (16) 

— OO 

where gi{e) = — Now, of course, the effective interaction in external channel is not 

reduced to sum of pairwise effective interactions with replacement e* —E — ti. Neverthe¬ 
less, this interaction includes three terms Wjk and is still suitable for Faddeev reduction. 
But now there are no pure pairwise forces (except n®^) in the effective Hamiltonian for 
the external three-body channel. 

Moreover, if even the external interaction Vi is disregarded at all, each term Wjk in 
the effective Hamiltonian (12) includes a dependence on the kinetic energy of the third 
particle, i.e., can be considered, generally speaking, as a three-body force. This dependence 
on the third particle momentum reduces the strength of the effective interaction between 
two other particles due to a specific energy dependence of the coupling constants (see 
below). Therefore, one can say that there are no pure two-body forces in the three- 
body system in such an approach, with the exception of that part of interaction which is 
included in (for NN system it is just the peripheral part of meson exchange). 

3. DRESSED-BAG MODEL FOR NN FORCES 

Here, we give a brief description of the two-component DBM for the NN interaction. 
The detailed description has been presented in our previous papers [8, 24]. (The effective- 
field theory description for the dybarion model of nuclear force has been also developed 
recently [38].) The main assumptions of the DBM are following: 

(i) interacting nucleons form at small and intermediate distances {t^n ~ 1 fm) a 
compound state - the dibaryon or six-quark bag dressed with tt, a, and p-fields; 

(ii) the coupling of the external NN channel with this state gives the basic attractive 
force between nucleons at intermediate and small distances, the a-dressed bag giving 
the main contribution. 
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Thus, nucleon-nucleon system can be found in two different phase states (channels): 
the NN phase and the dressed 6 q bag phase. In the NN (external) channel the system is 
described as two nncleons interacting via OBE; in the internal 6 q + a channel the system 
is treated as a 6 g bag snrronnded by the strong scalar-isoscalar a-field (a “dressed” bag)®). 
The external two-nncleon Hamiltonian inclndes the peripheral part of one-pion- and two- 
pion-exchange (OPE and TPE, respectively) interaction and Coulomb interaction: 

^ —t + |u + V I (^ith soft cutoff) ^ 

In the simplest version of DBM we used a pole approximation for the dressed-bag 
(internal) resolvent 5 :“: 

(-) 

where |q;) is the 6 g part of the wavefunction for the dressed bag and |k) represents the 
plane wave of the a-meson propagation. Here, Eai}^) is the total energy of the dressed 
bag: 

Ea{\E) =ma + e^{k), (17) 

where 

e„{k) = k‘^/2ma + uj^{k) c:^ma + k‘^/2ma, rfi^ = ^ ( 18 ) 

rua + rua 

uj„{k) = \/m 1 + is relativistic energy of a-meson, m^j and rria are masses of a-meson 
and 6 g bag, respectively. 

The effective interaction w{E) resnlted from the conpling of the external NN channel 
to the intermediate dressed-bag state is illustrated by the graph in Pig. 1. 



Figure 1. Effective NN interaction induced due to the production of an intermediate dressed bag. 


To derive the effective interaction w for NN channel in such an approximation, the 
knowledge of full internal Hamiltonian h™ of the dressed bag, as well as the full transition 

Full description of the NN interaction at energies £1 ~ 1 GeV still requires other fields in the bag, 
such as 27r, p, and w. 
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operator is not necessary. We need only to know how the transition operator acts 

on those dressed-bag states, which are included into the resolvent (16): k). The 

calculation of this quantity within a microscopical quark-meson model results in a sum 
of factorized terms [24]: 

= (19) 

L 

where 99 ;^^ G is the NN transition form factor and i?£(k) is the vertex function 
dependent on the a-meson momentum. 

Here we should elucidate our notation in respect to the quantum numbers of angular 
momenta. In general, the 6 g-state index a includes all the quantum numbers of the 
dressed bag, i.e., a = {J, M, S', T, Lf,, Lg-}, where Lj,, S', T, J, and M are the orbital 
angular momentum of the 6 q bag, its spin, isospin, total angular momentum, and its 
projection on the axis, respectively, and is the orbital angular momentum of the 
a meson. However, in the present version of the DBM, the s-wave state of the 6 q bag 
with the configuration only is taken into account, so that Lb = 0, J = S, and thus the 
isospin of the bag is uniquely determined by its spin. The states of the dressed bag with 
0 should lie higher than those with = 0. For this reason, the former states are 
not included in the present version of the model. Therefore, the state index a is specified 
here by the total angular momentum of the bag J and (if necessary) by its projection 
M: 

Thus, the effective interaction in the NN channel w{E) = can be 

written as a sum of separable terms in each partial wave: 


w{E) = wfi,(r,r',E), 

J,L,L' 


( 20 ) 


with 

r') = 5^ '(r) ^if'*(r'). (21) 

M 

The energy-dependent coupling constants A;^^/(i?) appearing in Eq. (21) are directly 
calculated from the loop diagram shown in Fig. 1; i.e., they are expressed in terms of 
the loop integral of the product of two transition vertices B and the convolution of two 
propagators for the meson and quark bag with respect to the momentum k: 



Bi(k)Bi;(k) 

E - E„(k) 


( 22 ) 
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The vertex form factors and the potential form factors G have been 

calculated in the microscopic quark-meson model [8, 24], 

When the iViV-channel wavefunction T™ is obtained by solving the Schrodinger equa¬ 
tion with the effective Hamiltonian the internal {6qN) component of the wave- 

function is found from Eq. (4): 


4<:;V(J5) = |a''">E-FV 


S^(k) 




E - E„(k) 


where the underlined part can be interpreted as the mesonic part of the dressed-bag 
wavefunction. 

The weight of the internal dressed-bag component of bound-state wavefunction (with 
given value J) is proportional to the norm of 


As one can see from the comparison between Eqs. (22) and (24), the integral in 
Eq.(24) is equal to the energy derivative (with opposite sign) of the coupling constant 

>~iu(E)- 


-'ll' — 


and thus we get an interesting result: 




dA(E) 
dE ’ 


i.e., the weight of internal QqN state is proportional to the energy derivative of the 
coupling constant of effective NN interaction. In other words, the stronger the energy 
dependence of the interaction in NN channel, the larger the weight of channel correspond¬ 
ing to non-nucleonic degrees of freedom. This result is in full agreement with well-known 
hypothesis: energy dependence of interaction is a sequence of underlying inner structure 
of interacting particles. 

The total wavefunction of the bound state T must be normalized to unity. Assum¬ 
ing that the external (nucleonic) part of the wavefunction found from the effective 
Schrodinger equation has the standard normalization ||T®""|| = 1, one obtains that the 
weight of internal, i.e., the dressed-bag component is equal to: 


Pi. = 




(26) 
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Thus, the NN interaction in DBM approach is a sum of peripheral terms and 

^TPE) representing OPE and TPE with soft cutoff parameter and an effective 

interaction w{E) (see Eqs. (20), (21)), which is expressed (in a single-pole approximation) 
as an one-term separable potential with the energy-dependent conpling constants (22). 
The potential form factors are taken as the conventional harmonic oscillator 

wavefunctions \2S) and \2Dy\ Therefore, the total NN potential in DBM model can be 
represented as: 

= pOPE -f -f + w{E) + XT, (26) 

where P = |(^o)(v^o| is the projector onto IQS') harmonic oscillator fnnction and constant 
A shonld be taken snfficiently large. 

Table 1. Deuteron properties in DBM and other current NN models 


Model 

Ed, MeV 

Pd, % 

rm, fm 

Qd, fm^ 

IXd, n.m. 

Asjfm 

V{D/S) 

RSC 

2.22461 

6.47 

1.957 

0.2796 

0.8429 

0.8776 

0.0262 

Moscow 99 

2.22452 

5.52 

1.966 

0.2722 

0.8483 

0.8844 

0.0255 

Bonn 2001 

2.224575 

4.85 

1.966 

0.270 

0.8521 

0.8846 

0.0256 

DBM (I) Peg = 3.66% 

2.22454 

5.22 

1.9715 

0.2754 

0.8548 

0.8864 

0.02588 

DBM (II) Peg = 2.5% 

2.22459 

5.31 

1.970 

0.2768 

0.8538 

0.8866 

0.0263 

Experiment 

2.224575 

- 

1.971 

0.2859 

0.8574 

0.8846 

0.0263'") 


'^^An average value of the asymptotic mixing parameter rj over the results of a few most accurate 
experiments is presented here (see [40-43]). 


The model described above gives a very good description for singlet ^S'o, the triplet 
— ^Di phase shifts, and mixing parameter ei in the energy region from zero np to 
1 GeV [24]. The denteron observables obtained in this model withont any additional or 
free parameter are presented in Table 1 in comparison with some other NN models and 
experimental values. The quality of agreement with experimental data for the NN phase 
shifts and deuteron static properties found with the presented force model, in general, 
is higher than those for the modern NN potential model snch as Bonn, Argonne, etc., 
especially for the asymptotic mixing parameter rj and the denteron qnadrnpole moment. 

It was first suggested [39] long ago and then confirmed in detailed 6q microscopic calculations [15] that 
the 6q wavefunction in NN channel corresponds just to excited 6 q-bag components | [42]TS'T), 

while the ground state |s®[ 6 ]) describes the wavefunction in the bag channel. 
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The weight of the internal (dressed-bag) component in the deuteron is varied from 2.5 to 
3.6% in different versions of the model [8, 24]. 

4. THREE-NUCLEON SYSTEM WITH DBM INTERACTION 

Eor description of the three-body system with the DBM interaction the momentnm 
representation is more appropriate. We will employ the same notation for fnnctions both 
in the coordinate and momentnm representations. The following notations for coordinates 
and momenta are employed: rj(pj) is relative coordinate (momentum) of pair (jk), while 
Pj (qj) is Jacobi coordinate (momentum) of *th particle relatively to the center of mass 
for the pair (jk), and k is usually a momentum of a meson. 


4.1. Effective interaction due to pairwise NN forces 


One obtains an effective Hamiltonian for the external 3N channel according to a general 
recipe for transition from two- to three-particle system: 

= T-h5^{nr-f-W(U)}, (27) 

i 

where each of three effective potentials takes the form: 

W(U) = J(q, - q')m,(U - gf/2m), (28) 


and fh = misima/{rriM -|- rria) is a reduced mass of nucleon and 6g bag. In the pole 
approximation, this effective interaction reduces to a sum of two-body separable potentials 
with the coupling constants depending on the total three-body energy E and the third- 
particle momentum gp 


= <^(q* 












(29) 


When using such an effective interaction, one must also include an additional 3BE due 
to the meson-exchange interaction between the dressed bag and the third nucleon (see 
the next subsection). The pattern of different interactions arising in the 3N system in 
such a way is illustrated in Eig. 2. 

In the single-pole approximation, the internal (dressed-bag) components of the total 
wavefunction are expressed in terms of the nucleonic component 'k®""(pj, q*) as 


vl/“(k,q,;U) = 




\a 


JiMi 




E — En — 


2m 


(30) 
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Figure 2. Different interactions in the iN system for one of three possible combinations (1 + 23) of 
three nucleons: the peripheral two-nucleon interaction vf^ is due to OPE -1- TPE, the effective two-body 
interaction Wi{E) is induced by the production of dressed 6q bag and meson-exchange 3BF . 


where are the overlap integrals of the external 3iV component and the potential 

form factors 



These overlap fnnctions depend on the momentnm (or coordinate), spin, and isospin 
of the third nncleon. For brevity, the spin-isospin parts of the overlap fnnctions and 
corresponding quantum numbers are omitted unless they are needed. In Eqs. (29)- 
(31) and below, we keep the index i in the quantum numbers Li and Ji in order to 
distinguish the orbital and total angular momenta attributed to the 2N form factors from 
the respective angular momenta J and L of the whole 3iV system. 

It should be noted that the angular part of the function (q*) in Eq. (31) is not 
equal to (<?). This part includes also other angular orbital momenta due to coupling 

of the angular momenta and spins of the dressed bag and those for the third nucleon. In 
the next section we consider the spin-angular and isospin parts of the overlaps functions 
X^(^*(qi) in more detail. 

The norm of each 6 qN component for the 3N bound state is determined by sum of 
the integrals: 




E 

JiMi 


a 


JiMi 





E — Erv — 


—V 

2m J 


■dk 




dq*. (32) 


The internal loop integral with respect to k in Eq. (32) (in braces) can be replaced by the 
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energy derivative of X'l: 


il£(k)il^Uk) 


E-E^ 


—V 

2m J 


dk 


J_\Ji 

'dE 



(33) 


Thus, the weight of the 6 qN component in the 3N system is determined by the same en¬ 
ergy dependence of the coupling constants as the contribution of the 6g component 

in the NN system but at a shifted energy. 

With using Eq. (33), the norm of QqN component can be rewritten eventually as 

JiMi LiL[ 



Due to explicit presence of the meson variables in our approach, it is generally impossi¬ 
ble to define the wavefunction describing relative motion of the third nucleon ^ipiq) in 
the 6 qN channel. However, by integrating T™(k, q) over the meson momentum k, one 
can obtain an average momentum distribution of the third nucleon in the 6 qN channel 
(i.e., those weighted with the cr-meson momentum distribution). Based on Eq. (33), we 
can attribute the meaning of the third nucleon wavefunction in the 6 qN channel to the 
quantity 



With this “quasi-wavefunction”, one can calculate the mean value of any operator depend¬ 
ing on the momentum (or coordinate) of the third nucleon. We note that the derivative 
—d\/dE is always positive. 


4.2. Three-body forces in the DBM 





Figure 3. The graphs corresponding to three new types of three-body force. 

In this study, we employ the effective interaction (29) and take into account the 
interaction between the dressed bag and the third nucleon as an additional 3BE. We 
consider here three types of 3BE: one-meson exchange (tt and a) between the dressed bag 
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and the third nucleon (see Figs. 3a and 3b) and the exchange by two a-mesons, where the 
third-nucleon propagator breaks the a-loop of the two-body force - 2(T-process (Fig. 3 c). 

All these forces can be represented in the effective Hamiltonian for external 3N channel 
as some integral operators with factorized kernels: 




E 

JM,J'M',L,L' 




3BF 




J'M' 


P* 


(36) 


Therefore, matrix elements for 3BF include only the overlap functions, and thus the 
contribution of 3BF is proportional to the weight of the internal 6 qN component in the 
total 3N wavefunction. To our knowledge, the first calculation of the 3BF contribution 
induced by OPE between the 6 q bag and the third nucleon was done by Fasano and 
Lee [44] in the hybrid QCB model using perturbation theory. They used the model where 
the weight of the Qq component in a deuteron was ca. 1.7%, and thus they obtained a very 
small value of-0.041 MeV for the 3BF OPE contribution to the 3N binding energy. Our 
results for the OPE 3BE agree with the results obtained by Easano and Lee (see Table 2 
in Section 7), because the OPE contribution to 3BE is proportional to the weight of the 6g 
component, and in our case, it should be at least twice as compared to their calculation. 
However, we found that a much larger contribution comes from scalar cr-meson exchanges: 
one-sigma exchange (OSE) and two-sigma exchange (TSE). We emphasize that, due to 
(proposed) restoration of chiral symmetry in our approach, the a-meson mass becomes 
ca. 400 MeV, and thus the effective radius of the a-exchange interaction is not so small 
as that in conventional OBE models. Therefore, we cannot use the perturbation theory 
anymore to estimate the 3BE contribution and have to do the full calculation including 
3BE in the total three-body Hamiltonian. 


4-2.1. One-meson exchange between the dressed bag and third nucleon 


Eor the one-meson exchange (OME) term, the three-body interaction takes 


the form: 


LM 


H^;Uq*,q'd^)= dk 


5i(k) 


E - Ea- qf/^m 


q') 


B^iy) 


E - Ea- qi/^m 


(37) 


Therefore, the matrix element for OME can be expressed in terms of the internal “bag” 
components 

(T""|OME|T'^") = (38) 
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The integral over a-meson momentum k (37) can be shown to be reduced to a difference 
of the values for constant \{E — q‘^/2m), so that the vertex functions B{k) can be excluded 
from formulas for OME 3BF matrix elements. The details of calculations for such matrix 
elements are given in the Appendix. 

4 -2.2. 2a-process 



Figure 4. The graph illustrating three-body scalar force due to two-sigma exchange (2(7-process). 

The 2(T-process (TSE) shown in Eig. 4 also contributes significantly to 3BE. This 3N 
interaction seems less important than the OSE force, because this interaction imposes a 
specific kinematic restriction on the 3N configuration®^. 

The operator of the TSE interaction includes explicitly the vertex functions for the 
transitions (NN 6q + <j) so that these vertices cannot be excluded similarly to the 

case of OME. Therefore, we have to choose some form for these functions. It is naturally 
to require that these vertices should be the same as those assumed in two-body DBM; 
i.e., they can be normalized by means of the coupling constants X{E), which, in turn, are 
chosen in the two-nucleon sector to accurately describe NN phase shifts and deuteron 
properties (see below Eq. (41) for vertex normalization). We use the Gaussian form factor 
for these vertices: 


^/2u;^{ky 

where k is the meson momentum and the parameter b is taken from the microscopical 


It follows from the intuitive picture of this interaction that this force can be large only if the momentum 
of the third nucleon is almost opposite to the momentum of the emitted a meson. Thus, a specific 
3N kinematic configuration is required when two nucleons approach close to each other to form a bag, 
while the third nucleon has a specific space localization and momentum. 
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quark model [24]: 


b 

Then, the vertex constants 



- —h^ 

- 24^0, 


hn = 0.5 fm. 


shonld be fonnd from the eqnation: 

= >^iAE), 




{E-ma-ea{k)) ■ 2 uJa{k) 


(40) 


(41) 


where are the conpling constants employed in the constrnction of the DBM in 

the 2N sector and are fixed by NN phase shifts. For the crNN vertices, we take also the 
Gaussian form factor: with = h^/Q. 

Then, the box diagram in Fig. 4 can be expressed in terms of the integral over the 
momentum qo of the third nucleon in the intermediate state: 


TSEti/JJ' 


W" (q, q'; B) = Sjj'Ann Bi‘- x 


X- 


dqo 




g-(a^+f'^){qo-q')^ 


(27r)3 J + (qo-q)^ E-m„-gg/2m ml + (qo-q')^’ 

Thus, the matrix element for the total contribution of TSE takes eventually the form 


(42) 


(TSE) = 3 ^ f 

J Ti /T T T I 'J 


x£^^(q) 


TSEtt/J iJi I 

uv} 


q.q'i-B)Xi; 


(43) 




After the partial wave decomposition, these six-dimensional integrals can be rednced to 
two-fold integrals, which are compnted nnmerically by means of the appropriate Ganssian 
quadratures. 

We should emphasize here that both two-nucleon force induced by the DBM and two 
parts of 3BF contribution in our approach, i.e., OSE and TSE, are all taken with unified 
conpling constants and nnified form factors in Eqs. (37), (39)-(41), in a sharp contrast 
to the traditional meson-exchange models (see also the section 8). 


5. VARIATIONAL CALCULATIONS OE 3N SYSTEM WITH DBM 

INTERACTION 


The effective Schrodinger eqnation for the external 3N part of the total wavefnnction 
^tot(^) 4 ,ex(^) _ with Hamiltonian 


3 

H^°\E) = T + + Wi{E) + W^^^{E)} 

i=l 


( 44 ) 
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has been solved by variational method using antisymmetrized Gaussian basis [45]. Be¬ 
cause of the explicit energy dependence of the three-body total Hamiltonian, we used an 
iterational procedure in respect to the total energy E for solving this equation: 

^tot/^(n-l)'v^ex(n) _ ^(n)^ex(n)^ 


Such iterations can be shown to converge, if the energy derivative of effective interaction is 
negative (for our case, this condition is valid always). For our calculations, 5-7 iterations 
provide usually the accuracy of 5 decimal digits for 3iV binding energy. 

Construction of a 3iV variational basis. Here, we give the form of the basis 
functions used in this work and the corresponding notation for the quantum numbers. 
The wavefunction of the external 3iV channel, can be written in the antisymmetrized 
basis as a sum of the three terms: 


^ex 


\]/(2) ^(3) 

ex ' ex ' ex 5 


(45) 


where the label {i) enumerates one of three possible set of the Jacobi coordinates (rj,/>J. 
Every term in Eq. (45) takes the form 

= («) 

7 n 

The basis functions are constructed from Gaussian functions and corresponding spin- 
angular and isospin factors: 


exp{-a,„rj - 




(47) 


where the spin-aiigular tind isospin components of the basis fnnctions are 

given in Appendix and the composite label 7 = 7(7 = {Xili L Sjk S tjk} represents the 
respective set of the quantum numbers for the basis functions (47): Aj is the orbital angular 
momentum of the (jk) pair; U is the orbital angular momentum of the third nucleon [i) 
relatively to the center of mass for the (jk) pair; L is the total orbital angular momentum 
of the 3N system; Sjk and tjk are the spin and isospin of the (jk) pair, respectively; and 
S is the total spin of the system. We omit here the total angular momentum J = 1/2 and 
its 2 :-projection M, as well as the total isospin of the system T = 1/2 and its projection 
Tz (in this work, we neglect the very small contribution of the T = 312 component). 

The nonlinear parameters of the basis functions a^n and f3^n are chosen on the Cheby- 
shev grid, which provides the completeness of the basis and fast convergence of variational 
calculations [46]. As was demonstrated earlier [47], this few-body Gaussian basis is very 
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flexible and can represent quite complicated few-body correlations. Therefore, it leads 
to the accurate eigenvalues and eigenfunctions. The formulas for the matrix elements of 
the Hamiltonian (for local NN interactions) on antisymmetrized Gaussian basis are given 
in the paper [45]. The matrix elements of DBM interactions on this basis are given in 
Appendix. 

Wavefunction in the internal 6 qN channel. Having the 3N component T 37 V found 
in the above variational calculation, one can construct the inner 6 giV-channel wavefunc¬ 
tion which depends on the coordinate (or momentum) of the third nucleon and the 
(T-meson momentum and includes the bag wavefunction (see Eq. (30)). By integrating 
the modulus squared of this function with respect to the meson momentum and inner 
variables of the bag, one obtains the density distribution of the third nucleon relative to 
the 6 q bag in the 6 qN channel. This density can be used to calculate further all observ¬ 
ables, whose operators depend on the variables of the nucleons and the bag. However, 
it is much more convenient and easier to deal with the quasi-wavefunction of the third 
nucleon in the 6 qN channel, which has been introduced by Eq. (35). 

To calculate matrix elements of the 3BE Coulomb and OPE forces, one needs the spin- 
isospin part of 6 qN components of the total wavefunction. Here we give them explicitly. 
The potential form factors now include the spin-isospin part (Pi)^? with 

quantum numbers corresponding to the dressed bag: 


Ji^i^d^dz 

^LiSd 






^d ’ 


= Itjtt ■■ tdtd.)- 


(48) 


The full set of the quantum numbers labelling the form factors includes the total (J*) and 
orbital (Lj) angular momenta, related to the vertex form factor, and also the spin and 
isospin numbers Sd, td, and td^ related to the dressed bag. However, since the present 
version of the DBM involves the bag states with zero orbital angular momentum, we have 
Sd = Jii while the bag spin and isospin are supplementary to each other: td + Sd = 1 . 
Hence we will omit the quantum numbers Sd and td, where they are unnecessary. 

The total overlap function Xi.^*(i) = [Tsat) can be written (with its spin-isospin 

part), e.g., as 


X£^(q4) = 


liJ 




Jmj 
1 ' 
2 


Hi 


{tdtd^ ntzi \TTz) Tu 


(49) 


Here, J and M are the total angular momentum of the 3N system and its ^(-projection, 
T and are the total isospin of the 3N system and its ^(-projection, while k and J' 
are the orbital and total angular momenta of the third (Ah) nucleon, respectively, and 
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is isospinor corresponding to the third nucleon. In the present calculation for the 
ground states of and ^He (with J = 1/2), we have considered two lowest even partial 
wave components {S and D) in 3iV wavefunctions only. Therefore, U can take only two 
values: 0 or 2. Moreover, the total angular momentum of the third nucleon J is uniquely 
determined by value of k: J = 1/2 at /* = 0 and J = 3/2 at A* = 2. So, actually there 
is no summation over J in Eq.(49). 

It is easy to see that the three form factors 99 ^) used in the present work ( 99 Q, and 
i^\) determine five radial components of the overlap function and five respective 

components of the quasi-wavefunction for the 6qN channel. To specify these components 
it is sufficient to give three quantum numbers, e.g., Sd, h and Lj, and we will use notation 
for these radial components: 

'^'oo,o: {Ji = Sd = td = 1, Li = 0, b = 0, = i), 

{Ji = Sd = l,td = 0, Li = 0, b = 0, J = i), 

(Li = Sd = l,td = 0, Li = 0,k = 2,J = |), 

^ 10 , 2 : {J^ = Sd = l,td = 0 , L, = 2,k = {),J = 1 ), 

2 : (L, = Sd = l,td = 0, L, = 2,k = 2,J = I). 

At last, we give a formula for the total quasi-wavefunction in internal channel (i), sep¬ 
arating out explicitly its spin-angular and isospin parts, which include the spin-isospin 
part of the bag wavefunction: 


=E r ■■ JM) m ■■ TT,)- 

\ Li 


( 60 ) 


The explicit dependence of this function on the isospin projection is important for 
calculation of the Coulomb matrix elements and r.m.s. charge radius. 

The interaction matrix elements include the overlap integrals of the potential form 
factors with the basis functions „ = <h 7 ^n + *h 7 ^n + *h 7 ^n, where all five above components 
of the overlap function enter the matrix elements independently (certainly, some of the 
matrix elements can vanish). The explicit formulas for the above overlap functions and 
detailed formulas for the matrix elements of all DBM interactions are given in Appendix. 
When calculating both the normalization of the internal components and observables, the 
QqN components distinguishing only by their radial parts can be summed. Thus, only 
three different components of the 6qN quasi-wavefunction remain: the one L-wave singlet 
{Sd = 0): 

\T/in = 

^00 — ^ 00 , 0 ) 
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and the two triplet ones {Sd = 1): 

^in _ ^in I ^in 
^10 ~ ^ 10,0 ^ 10 , 2 ! 

^in _ ^in I ^in 
^12 — ^ 12,0 ^ 12 , 2 - 

The total weight of each of three 6qN components is eqnal to 

d„'’ = Il'I'ioll" + ll't'.oll' + lll'ijf; i = 1.2,3. 


(61) 


(52) 


Now, let ns introdnce the relative weights of the individnal 6qN components: 


pm _ 
^SO — 


I Ip in I 
I ^ool 


P 


(*) 


pm _ 

^Sl — 


l^lol 


p 


(d 


pm _ 

— 




p 


(d 


(53) 


After renormalization of the fnll fonr-component wavefnnction, the total weight of all 
internal components is eqnal to 


P. = 

in 


3P 


(d 


1 + 3P 


(d 


(54) 


(here, we assnme that the 3iV component of the total wavefnnction, obtained from 
the variational calculation, is normalized to unity), while the total weight of the 3iV 
component Tsat is equal to 


P = 

p-v 


1 + 3K 


(d 


= 1 - P. 


(55) 


It is also interesting to find the total weight of the D wave with allowance for non- 
nucleonic components: 


Pd = P|r(l - Pm) + P'SPm. (56) 

Numerical values of all above probabilities for internal and external components are 
given below in Table 2. The total weight of all QqN components P^qN = Pm in the 3iV 
system turns out to be rather large and approaches or even exceeds 10%. Furthermore, 
taking into account the short-range character of these components, the more hard nucleon 
momentum distribution (closely associated with the first property) for these components, 
and very strong scalar three-body interaction in the internal QqN channels, one can con¬ 
clude that these non-nucleonic components are extremely important for the properties of 
nuclear systems. 


6. COULOMB EFFECTS IN ^He 

In this section we will demonstrate that the DBM approach leads to some new features 
related to the Coulomb effects in nuclei, and in particular in ^He. Eirst of all, the addi¬ 
tional Coulomb force arises because the 6g bag and rest nucleon can have electric charges. 



27 


We have found that this new Coulomb three-body force is responsible for a significant 
part of the total ^He Coulomb energy (this three-body Coulomb force has been missed 
fnlly in previons 3N calcnlations within hybrid 6qN models [48]). 

The second featnre of the interaction model nsed here is the absence of the local NN 
short-range repnlsive core. The role of this core is played by the condition of orthogonal¬ 
ity to the confined 6g states forbidden in the external NN channel. This orthogonality 
requirement imposed on the relative-motion NN wavefunction is responsible for the ap¬ 
pearance of some inner nodes and respective short-range loops in this wavefunction. These 
short-range nodes and loops lead to nnmerons effects and general conseqnences for the 
nnclear strnctnre. One of these conseqnences is a rather strong overestimation of the 
Conlomb contribntion when nsing the interaction between point-like nncleons. Thns, it 
is necessary to take into account the finite radius of the nucleon charge distribution^^ 

At last, in order to obtain the accurate Coulomb displacement energy AEq = 
— Esi^He), one should take into consideration the effects associated with the 
small mass difference between the proton and nentron. It is well known [49] that the 
above mass difference makes rather small contribntion to the difference between ^He and 
binding energies. Therefore, it was taken nsnally into acconnt in a pertnrbation ap¬ 
proach. However, since the average kinetic energy in our case is twice the kinetic energy 
in conventional force models, this correction is expected also to be much larger in our 
case. Hence, we present here the estimation for such a correction term without usage of 
the perturbation theory. 


6.1. “Smeared” Coulomb interaction 


The Ganssian charge distribntion p(r), that corresponds to the r.m.s. charge radins r, 
and is normalized to the total charge 2 :: 47r f pr'^dr = can be written as 


P(r)=z[^) 


a-^ = V. 

3 


(57) 


The Conlomb potential for the interaction between snch a charge distribntion p{r) and a 
point-like charged particle has the well-known form 


V{R) = 


dr p(r) 
|R — r| 


^erf {R^/a). 
K 


9 ) 


We remind in this point that the account of the finite radii of nucleons in the conventional approaches 


leads to fully negligible corrections to the Coulomb energy. 
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We have derived here a similar formula for the Coulomb interaction between two Gaussian 
distributions with different widths ai and 0:2 and r.m.s. radii and Tcj, respectively: 

V{R]ai,a 2 ) = ^^erf(i?yd), a = , or a~^ = +rl^)- (58) 

it (y.\ + 0^2 o 

In our calculations, we used the following charge radii for the nucleon and dibaryon: 

{rc)p = 0.87 fm, 

(rc)6g = 0.6 fmi°). 

These values lead to the “smeared” Coulomb interactions in the NN and QqN channels: 

V) = ^erf(rv/aAr7v), = 1.005 fm, 

= 0.863 fm. (59) 

6.2. Matrix elements of the three-body Coulomb force 


The Conlomb interaction between the charged bag and the third nncleon in the 
6qN channel is determined by the three-particle operator with the separable kernel (see 
Eq. (36)): 

c-'W«(p„p';q„q') = 


JiMiLiL'. 




Pi 


l+r. 


(b 




JiMi 


( 60 ) 


where (l-|-r 3*^)/2 is the operator of the ith nucleon charge and 1 + is operator of the 
bag charge. It is evident that the matrix element of the operator (60) can be expressed in 
terms of the integrals of the prodnct of the overlap fnnctions of form factors 

and three-body basis fnnctions. The method for calcnlation of snch Conlomb integrals is 
given in Appendix. 


This value is simply the r.m.s. charge radius of the 6q bag with the parameters given in [24]. The 
neutral a field of the bag changes this value only slightly. The evident difference between the charge 
radii of the nucleon and dibaryon can be well understood as follows: the charge radius of the 3q core 
of the nucleon is taken usually as ~ 0.5-0.55 fm, while remaining 0.3 fm is assumed to come from 
the charge distribution of the 7r+ cloud surrounding the 3q core in the proton. In contrast, the meson 
cloud of the dibaryon in our approach is mainly due to the neutral scalar-isoscalar a meson, so that 
the dibaryon charge distribution is characterized by the charge radius of the bare 6q core only. 
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7. RESULTS OF CALCULATIONS 

Here, we present the results of the 3N bound-state calculations based on two variants 
of the DBM. 

(I) In the first version developed in [24], the dressed-bag propagator includes three 
loops (two loops are with pions and one loop is with a meson), two of them are of 
the type shown in Fig. 2 of [24], in which each loop was calculated within the ^Pq 
model for quark-meson interaction. The third loop consists of two such vertices 
and a convolution of the a-meson and 6g-bag propagators [24]. 

(II) In the second version, we replaced two above pionic loops with the effective Gaussian 
form factor B{k), which describes the direct NN 6q +a transition, i.e., the direct 

transition from the NN channel to the dressed-dibaryon channel. 

Both versions have been fitted to the NN phase shifts in low partial waves up to an 
energy 1 GeV with almost the same quality. Therefore, they can be considered on equal 
footing. However, version (H) has one important advantage. Here, the energy dependence 
arising from the convolution of the two propagators involved into the loop, i.e., the 
propagators of the a meson and bare dibaryon, describes (with no further correction) just 
the energy dependence of the effective strength of the NN potential which is 

thereby taken directly from the above loop integral. In contrast, in the first version of the 
model, two additional qqnna loops give a rather singular three-dimensional integral for 
A‘'^^(U), where the energy dependence at higher energies should be corrected by a linear 
term. 

7.1. Bound-state energies of^H and ^He and individual contributions to them 

The main difference between the results for both versions is that the energy dependence 
of \{E) for the second version is much weaker than that for the first one. In addition, this 
energy dependence leads to some decrease in the contribution of the 6qN component to all 
3N observables and thus to the respective increase of the two-body force contribution as 
compared to the three-body force one. Table 2 presents the calculation results for the two 
above versions for the following characteristics: the weights of the internal 6qN channels 
and D wave in the total 3N function, as well as the weight of the mixed-symmetry 
S' component (only for the 3N channel); the average individual contributions from the 
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Table 2. Results of the 3N calculations with two- and three-body forces for two variants of the DBM 


Model 

E, MeV 

Pd, % 

Ps', % 

PeqNiPin), % 

Contributions to H, MeV 

T 

T + 

y(3iV) 


DBM(I) g = 9.577^^') 

-8.482 

6.87 

0.67 

10.99 

112.8 

-1.33 

-7.15 

DBM(II) g = 8.673('") 

-8.481 

7.08 

0.68 

7.39 

112.4 

-3.79 

-4.69 

AV18 + UIXC^) 

-8.48 

9.3 

1.05 

- 

51.4 

-7.27 

-1.19 

^He 

DBM(I) 

-7.772 

6.85 

0.74 

10.80 

110.2 

-0.90 

-6.88 

DBM(II) 

-7.789 

7.06 

0.75 

7.26 

109.9 

-3.28 

-4.51 

AV18 + UIXC^) 

-7.76 

9.25 

1.24 

- 

50.6 

-6.54 

-1.17 


^^These values of aNN coupling constant in calculations have been chosen to reproduce the exact 
binding energy of nucleus. The calculations for ^He have been carried out without any free 
parameters. 

^’^The values are taken from [50]. 


kinetic energy T, two-body interactions pins the kinetic energy T, and three-body 

force dne to OSE and TSE to the total Hamiltonian expectation. 

Eor variant I of the model, we present also the resnlt calcnlated when both 3BE and 
the dependence of the effective two-body force on the momentum of the third nucleon 
are disregarded (the first line). The results in the second line of Table 2 are obtained 
including the dependence of pair forces, but disregarding 3BE. The percentages of the 
D-wave and the internal components given in Table 2 were obtained with incorporation 
of the three internal components; i.e., these values correspond to the normalization of the 
total (four-component) wavefunction of the system to unity. 

To compare the predictions of the new model with the respective results for the con¬ 
ventional NN potential models, Table 2 also presents the results of recent calculations 
with the Argonne potential AV18 and Urbanna-Illinois 3BE UIX [50]. 

7.2. The densities, r.m.s. radii and charge distributions in and ^He 
At first, we give definitions of the nucleon and charge distributions in multichannel 


system. 
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The external 3N channel. The proton (pp) and neutron (p„) densities in this channel 
are defined by the standard way [51]: 




r = 








i=l 


- |P*) 1 ± T-; 


(b 


\ = 


iV, 




il} 


- fPl) 1 ± T-; 


( 1 ) 




(61) 

where pi is Jacobi coordinate in the set (i) and N^pj is the number of protons (neutrons). 
Due to property 




E 

2=1 


1 ± 


(b 




l±n 


d/-) = |±T 3 = iV{.}, 


the above densities are normalized to unity, provided that the external wavefunction 
is also normalized to unity: 

J Pp>}{r)r'^dr = 1 . 

The matrix element is proportional to 2 :-projection of the total isospin 

T 3 , therefore, the nucleon densities can be separated into isoscalar (matter) density and 
isovector parts: 


psir) = pm = 
Pv{r) = ^ ^ 


J(r - fpi) 




^{r-^pi) ( 1 ) 

p2 ^3 




(62) 


(63) 


Both latter densities are also normalized to unity. Then the nucleon densities can be 
expressed in terms of isoscalar and isovector densities as: 


p;^(^He)=p“(^H) = -(3p, + p,), 


p-(^H)=p-(^He) = -(3p,-p„). 
R.m.s. radii of corresponding distributions are equal to: 


(64) 


(’■ = / PTs.v.r.n](r)r dr. 


( 66 ) 


The r.m.s. charge radius in the 3N sector is also defined conventionally: 

(bh>“ = {'X + K + (“> 

where = 0.7569 fm^ and R^ = —0.1161 fm^ are the squared charge radii of the proton 
and neutron, respectively. 

The various types of one-particle densities (isoscalar, isovector, proton, neutron) in 
external 3N channel for the and ^He ground states calculated in DBM(I) are shown 
in Fig. 5. 
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r, fm 


Figure 5. The isoscalar ps, isovector proton pp, and neutron p„ densities in external 3N channel 
for and ^He systems obtained with DBM (version I). 


Below we present also two-proton density for ^He, which is defined nsnally as [52]: 

. « /mex 1 - ri) 1 + rf 1 + rf I \ 


Pppir) = 6 (^vl/ex - y -. (67) 

This density is normalized to 2: f Ppp{r)r‘^dr = 2. (As there is only a single nncleon in 
6qN channel, we do not attach the index “ex” to this quantity.) The two-neutron density 
Ppp{r) for is defined similarly (with replacing 1 -|- 1 — in Eq. (67). In Fig. 6 

we show both these densities for DBM (I) and also the two-proton density for ^He found 
with Bonn NN potential [52]. 

The internal 6qN channels. Here we define a density (normalized to unity) of the 


pure nucleon distributions as 


/win pVi 

^ n J 


2_/min S{r - api)l±T^ 

pW \ ^ 2 


T- , 


where = (T™!'!'™) and the quantity 


/wui — _ 


,in 1±< 




has the meaning of the average number of protons (neutrons) in the one internal 6qN 
channel (note that N^ + N^ = 1, i.e., there is only one nucleon in each internal channel). 
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Figure 6. The two-proton density in ^He (solid line) and two-neutron density in (dashed line) 
calculated with DBM (version I) in comparison with two-proton density found with Bonn potential [52] 
(triangles). 


The number iVl“ , depends on ratio of norms of 6qN components with different values of 

^ n J 

isospin of the bag. Therefore, the separation of the 6giV-channel density into isoscalar 
and isovector parts has no meaning. 

These average numbers of nucleons in 6qN channel can be expressed through relative 
probabilities of the 6qN components with definite value of isospin t, which in our case 
are equal to: 


pin _ pin _ pin i pin 

0 — ^t=0 — ^S1 5 

pin — pin _ pin 
— ^t=l — -^ 50 ’ 


(70) 


where P^, and are determined by Eq. (53). Hence, P/dg + P^bi = 1. Then one 
can write down the average numbers of nucleons as 

= Af”(’He) = 

Af“(’He) = ]V“(’H) = + if;”,. (71) 

The nucleon densities (61) can be expressed by similar formula through components of 
the internal wavefunction with definite value of isospin t. 

The total densities of nucleon distributions. The total nucleon densities (normal¬ 
ized to unity) for whole 3N system with allowance for both the 3N and 6qN components 
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can be now defined as 


Pil} 


(1 - Pin)pfp.iV{P} + 

bn-f ^ b7T,J b7j,J 

Win) 


(72) 


where Pin = 3P™/(1 + 3P™) is the total weight of all three internal channels (Eq. (54)) 
and the denominator: 


Wn) = (1 - ( 73 ) 

is equal to the average number of protons (neutrons) in the whole multicomponent system. 
The densities of the total proton and neutron distributions and also external- and internal- 
channel distributions for ^He calculated for DBM(i) are presented in Fig. 7. 



Figure 7. The external internal p™, and total p(total) densities of proton- and neutron 
distributions in ^He found with DBM (version I). 

One can define also a (normalized) density of the matter (or mass) distribution in the 
QqN channel as: 

r'plP) = + ■)('• - (1 - a)p,)m,|4-;“), (74) 

where a = md/{jnd -|- ttiat), rriN is a nucleon mass, and rud is mass of the bag (dibaryon). 
Then the total matter density (normalized to unity) is equal to 

(1 - Pin)p^3mAr -h Pir^pZi^N + nid) 

{m) 


r^Pm(r) 


(75) 
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The r.m.s. radius of any distribution normalized to unity is defined by Eq. (65). The 
denominator in Eq. (75) determines the average mass of the whole system with taking 
into acconnt non-nncleonic channels: 

(m) = (1 - Pin)3mAr + PinirriN + md) = Stun + Pmirrid - 2mN) > Stun- (76) 

Table 3. Isospin structure of 6qN channel, average number of nucleons and average mass calculated 
with the ground state wavefunctions of and ^He in the DBM approach 



3 

H 

^He 


DBM(I) 

DBM(II) 

DBM(I) 

DBM(II) 

pin 

^t=l 

0.6004 

0.6005 

0.6044 

0.6044 

pin 

0.3996 

0.3995 

0.3956 

0.3956 

CO 

0.799 

0.799 

2.209 

2.209 

3iV)f 

2.201 

2.201 

0.791 

0.791 

w 

0.919 

0.945 

1.863 

1.908 

(Pn) 

1.861 

1.906 

0.920 

0.947 

{N) 

2.780 

2.852 

2.784 

2.855 

{m)/3m N 

1.015 

1.010 

1.014 

1.010 


In the Table 3 we present some characteristics of isospin structure for wavefunctions 
in the 6qN channel: the relative probabilities for the components with t = 0 and t = 1 
(i.e., and P^Ti), average nnmbers of protons and nentrons in all three 6qN compo¬ 
nents (3(iV{P})), and also the average nnmber of nncleons (N) and the average mass (m) 
(divided by valne) in the whole fonr-component 3N system. It shonld be noted that 
the average number of nucleons in our multicomponent model, (iV{P}), is always less than 
the numbers of nucleons in 3N channel just due to existence of the non-nucleonic compo¬ 
nents. Eor example for DBM (I), the average number of protons in is approximately 
eqnal to the average nnmber of nentrons in ^He, viz. (iVp)(^H) ^ and is eqnal 

0.92 while the average nnmber of nentrons in is approximately eqnal to the average 
number of protons in ^He, viz. (iV„)(^H) (iVp)(^He), and is equal 1.86. Hence the 

average number of nucleons found with the total multicomponent and ^He functions 
is also always less than 3: 

(N) = (N,) + (N^) = 3 - 2Pi, < 3. 

In our DBM, (N) is equal 2.78 and 2.85 for versions I and II, respectively. 
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The charge distributions. The charge distribution for the point-like particles in the 
6qN channel can be written as the charge density of a system consisting of a point-like 
nncleon and a point-like bag: 





6{r - a pi) -k (5(r - (1 - a)pi){l + 4 ) 


where l-l-ts is operator of the bag charge. The total charge radius in QqN channel includes 
the r.m.s. radius of this point-like distribution (?"^)ch-point) the nucleon charge radius (i?p 
or Rn) and the charge bag radius Rd, which depends on the bag isospin t and its projection 
ts- 



in _ 

ch 



in 

ch—point 





y^^tdaRditRs) 



(78) 


The last term in Eq. (78) inclndes the projectors T^ <3 onto the 6 g bag isospin state with 
definite valnes of isospin t and its projection ts and is eqnal to (for the and ^He states 
with total isospin T = 1/2): 


= - <1 RliO, 0)P£„ + 




pm 

^t=l 


(79) 

The own charge radius of the 6 g bag R^itRR has, in general, different values in different 
isospin states, (which is related to the different multiquark dynamics in the channels with 
different isospin valnes) bnt we snppose in this work that their difference can be ignored, 
viz. 


i?2(0, 0) = Rl{R 0) = Rl{l, l) = Rl = hi ho = 0.6 fm, 
and Rdi^, — 1 ) = 0 (which corresponds to a nn-bag). (80) 

With the above assumptions the 6 g bag contribution to the and ^He charge radius is 
reduced to: 

A(''d.)L« = § (fi" + ^-^Pi“) ■ ( 81 ) 

The r.m.s. charge radius of whole multicomponent system is defined as: 

rl = {l-RR{rl2 + nn{rltR 

In the Table 4 we give r.m.s. radii for all the above distributions in and ^He found 
in the impulse approximation, as well as the respective experimental values and results 
obtained for AV18(2iV) + UIX(3iV) forces. To demonstrate the separate contributions 
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Table 4. The total r.m.s. radii (in fm) for the proton (r^), neutron (r„), matter (r^), and charge 
(rch) distributions in DBM approach and their separate values for external and internal channels 


Model 

3H 

^He 

rp 

rn 

rm 

rch 

rp 

Tn 

rm 

rch 

DBM(I) 

3N 

1.625 

1.770 

1.723 

1.779 

1.805 

1.648 

1.754 

1.989 


6qN 

1.608 

1.823 

1.142 

1.188 

1.854 

1.618 

1.159 

1.412 


Total 

1.625 

1.773 

1.663 

1.724 

1.807 

1.647 

1.694 

1.935 

DBM(II) 

3N 

1.613 

1.761 

1.713 

1.769 

1.795 

1.636 

1.744 

1.980 


6qN 

1.573 

1.797 

1.124 

1.171 

1.829 

1.583 

1.141 

1.396 


Total 

1.613 

1.762 

1.672 

1.732 

1.796 

1.635 

1.703 

1.944 

AV18 -b 

1.59 

1.73 



1.76 

1.61 



Experiment 

1.60(^) 



1.755 

1 . 77 (b) 



1.95 


Taken from [53]. 

These “experimental” values are taken from [53]. They have been obtained by substruction of the own 
proton and neutron charge radii squared (0.743 and -0.116 fm^, respectively) from the experimental 
values of the charge radii squared. 

of the 3N and 6qN channels to these observables, we also present the valnes calcnlated 
separately with only nncleonic and 6qN parts of the total wavefnnction. It is seen from 
the Table 4 that both versions of onr model (viz. DBM(I) and DBM(II)) give qnite 
similar values for all the radii. The most interesting point here is the importance of 6qN 
component contributions. In fact, the contribution of the 6qN channel shifts all the radii, 
i.e., Teh and Vp in and ^He, predicted with pure nucleonic components in our approach, 
much closer to the respective experimental values. For example, the value rch = 1-822 fm 
calculated for with only the nucleonic part of the wavefnnction is essentially larger 
than the respective experimental value 1.755 fm. However, an admixture of a rather 
compact 6qN component (rch = 1-22 fm) immediately shifts the charge radius down 
to a value of 1.766 fm, which is very close to its experimental value. 

Thus, the dibaryon-nucleon component works in a right way also in this aspect. It 
is interesting to note that, in general, the predictions of our two-phase model are quite 
close to those of the conventional pure nucleonic AV18 + UIX model. This means that 
our multichannel model is effectively similar to a conventional purely nucleonic model (at 
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least for many static characteristics). However, this similarity will surely hold only for 
the characteristics that are sensitive mainly to low-momentum transfers, while the prop¬ 
erties and processes involving high-momentum transfers will be treated in two alternative 
approaches in completely different ways. 

7.3. Coulomb displacement energy and charge symmetry breaking effects 

The problem of accurate description of Coulomb effects in ^He in the current 3N 
approach of the Faddeev or variational type has attracted much attention for last three 
decades (see, e.g., [49, 54] and the references therein to the earlier works). It is interesting 
that the Coulomb puzzle in ^He, being related to the long-range interactions, is treated 
in a different manner in our and conventional approaches. 

The AEq problem dates back to the first accurate 3N calculations performed on the 
basis of the Faddeev equations with realistic NN interactions in the mid-1970s [55]. 
These pioneer calculations first exhibited a hardly removable difference of ca. 120 keV 
between the theoretical prediction for AE^ ~ 640 keV and the respective experimental 
value AE^^ ~ 760 keV. In subsequent 30 years, numerous accurate 3N calculations 
have been performed over the world using many approaches, but this puzzle was still 
generally unsolved. The most plausible quantitative explanation (but yet not free of 
serious questions) for the puzzle has been recently suggested by Nogga et al. [49]. They 
have observed that the difference in the singlet scattering lengths of pp (nuclear part) 
and nn systems (originating from the effects of charge symmetry breaking (CSB)) can 
increase the energy difference between and ^He binding energies and thus contribute 
to AEq. 

Our results obtained in this work with DBM give an alternative explanation of the 
AEq puzzle and other Coulomb effects in ^He without any free parameter. The Coulomb 
displacement energies AEq^ together with the individual contributions to the Ai7c-value, 
are presented in Table 5. 

We emphasize three important points, where our results differ from those for conven¬ 
tional models. 

(i) First, we found a serious difference between conventional and our approaches in the 
short-range behavior of wavefunctions even in the nucleonic channel. Conventional 
3N wavefunctions are strongly suppressed along all three interparticle coordinates 
Tij due to the short-range local repulsive core, while our wavefunctions (in the 3N 
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Table 5. Contribution of various terms (in keV) of the Coulomb interaction to the ^H-^He mass 
difference /S.Ec 


Contribution 

DBM(I) 

DBM(II) 

AV18 + UIX 

Point Coulomb 2>N only 

598 

630 

677 

Point Coulomb 2>N + 6qN 

840 

782 

- 

Smeared Coulomb 3N only 

547 

579 

648 

Smeared Coulomb 3N -|- 6qN 

710 

692 

- 

np mass difference 

46 

45 

14 

Nuclear CSB (see Table 6) 

0 

0 

-k65 

Magnetic moments and spin-orbit®"^ 

17 

17 

17 

Total 

773 

754 

754 


'^^Here we use the value for this correction from [49]. 

channel) have stationary nodes and short-range loops along all and the third 
Jacobi coordinates pk. Snch a node along the p coordinate is seen also in the 
QqN relative-motion wavefnnction. This very pecnliar short-range behavior of onr 
wavefnnctions leads to a strong enhancement of the high-momentnm components 
of nuclear wavefnnctions, which is required by various modern experiments. On the 
other hand, these short-range radial loops lead to significant errors, when using the 
Coulomb interaction between point-like particles within our approach. Hence, we 
must take into account the finite radii of charge distributions in the proton and 6g 
bag. Otherwise, all Coulomb energies will be overestimated. 

(ii) Another important effect following from our calculations is a quite significant con¬ 
tribution of the internal QqN component to /S.Eq. In fact, just this interaction, 
which is completely missing in conventional nuclear force models, makes the main 
contribution (163 and 113 keV for versions I and II correspondingly) to filling the 
gap in AT’c between conventional 3iV calculations and experiment if the CSB effects 
are of little significance in AEq- 

The large magnitude of this three-body Coulomb force contribution in our mod¬ 
els can be explained by two factors: first, a rather short average distance 
between the 6g bag and the third nucleon (which enhances the Coulomb interaction 
in the QqN channels) and, second, a relatively large weight of the QqN components, 
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where the 6q bag has the charge +1 (i.e., it is formed from np pair). This specific 
Coulomb repulsion in the 6qN channel should appear also in all other nuclei, where 
the total weight of such components is abont 10% and higher. Therefore, it shonld 
strongly contribnte to the Conlomb displacement energies over the entire periodic 
table and conld somehow explain the long-term Nolen-Schiffer paradox [56] in this 
way. 


(iii) The third specific effect that has been fonnd in this stndy and contribntes to the 
qnantitative explanation of AEq is a strong increase in the average kinetic energy 
(T) of the system. This increase in (T) has been already discovered in the first 
early 3N calculations with the Moscow NN potential model [57] and results in a 
similar nodal wavefunction behavior along all interparticle coordinates but without 
any non-nncleonic component. 

The increase in (T) leads to the proportional increase in the np mass difference 
correction to AEq. Since the average kinetic energy in onr case is twice the kinetic 
energy in conventional force models, this correction is expected to be also mnch 
larger in our case. Hence, we evaluate such a correction term in the following way 
(without usage of the perturbation theory). In the conventional isospin formalism, 
one can assume that the and ^He nuclei consist of the equal-mass nucleons: 

_ mp +TUn 


SO that nip = m + Am/2, mn = m — Am/2, where Am = mp — mn- The simplest 
way to inclnde the correction dne to the mass difference Am is to assnme that all 
particles in have the average mass 


mu 


2m„ -|- m. 


- - = m - Am, 

3 6 


while in ^He they have the different average mass 


mue = 


2m„ -|- m, 


= m -\ — Am. 
6 


In spite of smallness of the parameter Am/m, the pertnrbation theory in respect of 
this parameter does not work. So we nsed the average mass rhu in calcnlation of 
and friue in calculation of ^He. The contribution of this np mass difference to 
the AEq value is given in the fifth row of Table 5. As is seen from the table, this 
correction is not very small in our case and contributes to AEq quite significantly. 
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Many other effects attributed to increasing the average kinetic energy of the sys¬ 
tem will arise in our approach, e.g., numerous effects associated with the enhanced 
Fermi motion of nucleons in nuclei. 

Charge symmetry breaking effects in DBM. As was noted above, the best ex¬ 
planation for the Ai?c value in the framework of conventional force models published up 
to date [49] is based on the introduction of some CSB effect, i.e., the difference between 
nn and pp strong interactions. At present, two alternative values of the nn scattering 
length are assumed: 

®nn = —18.7 fm and a® = —16.3 fm. (82) 

The first value has been extracted from the previous analysis of experiments 
d{n~,'y)nn [58] (see also [59] and references therein) and is used in all current NN poten¬ 
tial models, while the second value in (82) has been derived from numerous three-body 
breakup experiments n + d ^ nnp done for the last three decades. In recent years, such 
breakup experiments are usually treated in the complete Faddeev formalism, which in- 
eludes most accurately both two-body and 3BF [60]. Thus, this aiin value is considered 
as a quite reliable one. However, the quantitative explanation for the AEc value in 
conventional force models uses just the first value of ann as an essential point of all the 
construction. At the same time, the use of the second value a„n(= —16.3 fm) (which is 
not less reliable than the first one) invalidates completely the above explanation! 

Therefore, in order to understand the situation more deeply and to determine the 
degree of sensitivity of our prediction for AEc to variation in ann, we made also 3N 
calculations with two possible values of a^n from Eq. (82). These calculations have been 
carried out with the effective values of the singlet-channel coupling constant corresponding 
to the VATgAT part of the NN force: 


-^3He( So) — -\p + -Anp, 


A3h(^>S'o) — —\nn + TjA^p. 


3' ' 3' 


(83) 

(84) 


In the above calculations, we employ the value \np = 328.9 MeV that provides the 
accurate description of the np phase shifts and the experimental value of the np 
scattering length = —23.74 fm [24]. Here, for pp-channel we use the value App = 
325.523 MeV fitted to the well-known experimental magnitude Opp = —8.72 fm and for 
nn-channel two \nn values corresponding to two available alternative values of the nn 
scattering length (82) have been tested. The calculation results are presented in Table 6. 
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Table 6. Contribution of charge symmetry breaking effects to the ^H-^He mass difference Ai?c 



AEc 

keV 

ann, fm 

DBM(I) 

DBM(II) 

-16.3 

-18 

-39 

-18.9 

+45 

+26 


As is seen in Tables 5 and 6, the DBM (version I) can precisely reproduce the Coulomb 
displacement energy AEq with the lower (in modulus) value ttnn = —16.3 fm, while this 
model overestimates AEc by 54 keV (=45 + 9 keV) with the larger (in modulus) value 
dnn = —18.9 fm. Thus, the DBM approach, in contrast to the conventional force models, 
prefers the lower (in modulus) possible value -16.3 fm of the nn scattering length, which 
has been extracted from very numerous 3N breakup experiments n + d ^ nnp [60]. 

Now, let us discuss shortly the magnitude of CSB effects in our model. The measure 
of CSB effects at low energies is used to consider the difference between a^n and so- 
called “pure nuclear” pp scattering length that is found from pp scattering data, when 
the Coulomb potential is disregarded. The model dependence of the latter quantity 
was actively discussed in the 1970s-1980s [61-63]. However, the majority of modern 
NN potentials fitted to the experimental value app = —8.72 fm results in the value 
ttpp = —17.3 fm, when the Coulomb interaction is discarded. It is just the value that is 
adopted now as an “empirical” value of the pp scattering length [64]. Thus, the difference 
between this value and is usually considered as the measure of CSB effects. However, 
our model (also fitted to the same experimental value app = —8.72 fm) gives a quite 
surprising result: 

a^(DBM) = -16.57 fm, (85) 

which differs significantly from the above conventional value (by 0.8 fm) due to the 
explicit energy dependence of the NN force in our approach. 

Thus, if the difference is still taken as the measure of CSB effects, the 

smallness of this difference obtained in our model testifies to a small magnitude of the 
CSB effects, which is remarkably smaller than the values derived from conventional OBE 
models for the NN force. 
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8. DISCUSSION 

The 3N results presented in the previous section differ significantly from the results 
found with any conventional model for NN and 3N forces (based on Ynkawa’s meson 
exchange mechanism) and also from the resnlts obtained in the framework of hybrid 
models [65], which inclnde the two-component representation of the NN wavefnnction 
4/ = Tjvtv + 'I'eg- It is convenient to discuss these differences in the following order. 

(i) We fonnd that the dependence of pair NN forces on the momentnm of the third 
particle in the 3N system is more prononnced in onr case than in other hybrid 
models [34, 48, 65, 66]: the 3N binding energy decreases by ca. 1.7 MeV, from 
5.83 to 4.14 MeV when one takes into account the g^-dependence (cf. the first 
and second rows in Table 2). From more general point of view, it means that, in 
onr approach, pairwise NN interactions (except Ynkawa OPE and TPE terms), 
being “embedded” into a many-body system, loose their two-particle character and 
become substantially many-body forces (i.e., depending on the momenta of other 
particles of the system). 

(ii) Dne to snch a strong g^ dependence (of “repnlsive” character), the 3N system cal¬ 
culated including only the pairwise forces turns out to be strongly underbound 
{E = —4.14 MeV). In other words, the “pairwise” NN forces (including their g^ 
dependence on the momenta of the third nucleon) give only about half the total 3N 
binding energy, leaving the second half for the 3BE contribntion. Therefore, the 
following qnestion is decisively important: can the 3BE (inevitably arising in onr 
approach) give the large missing contribntion to the 3N binding energy? Usefnlness 
of the developed model for the description of nuclear systems depends directly on 
the answer to this important question. It is appropriate here to remind that in the 
conventional 3BE models such as Urbana-Illinois or Tucson-Melbourne, the con¬ 
tribution of 3BE to the total 3N binding energy does not exceed 1 MeV; i.e., this 
contribntion can be considered as some correction (~ 15%), althongh it is significant 
for the precise description of the 3N system. 

(iii) Eortnnately, the contribntion of 3BE indnced by OSE and TSE enables one to 
fill this 4.3 MeV gap between the two-body force contribntion and experimental 
valne. In fact, inclnding both OSE and TSE contribntions to 3BE, taken with 
the same coupling constants and form factors as in the driving iViV-force model. 
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together with a quite reasonable value for the aNN coupling constant, q^nn = 8- 
10, one obtains the 3N binding energy that is very close to the experimental value 
(see rows 3, 4 and 6, 7 in Table 2). Thns, the presented force model leads to a 
very reasonable binding energy for the 3N system, however, with the mnrch larger 
(as compared to the traditional 3N force model) contribntion of 3BF. In fact, the 
unification of the basic 2N and 3N force parameters provides a strong support for 
the whole force model suggested here and is in a sharp contrast with all traditional 
force model based on t-channel exchange mechanism. We remind to reader that 
the 2N and 3N forces in conventional approaches (where the latter is indnced by 
an intermediate A-isobar prodnction) are taken with different cnt-off parameters 
valnes and A^tva in 2A and 3N sectors in order to explain the basic featnres 

of 3N nuclei and N + d scattering! Thus in the traditional approach, one has some 
serious inconsistency in parameter values for 2N and 3N sectors. 

(iv) The contribntions of the pairwise and different three-body forces to the total 3N 
binding energy for are given in the fifth and sixth rows of Table 2. From the 
results presented in this table, one can conclude that just the total 3BF contribution 
to the 3N binding energy dominates and, in fact, determines the whole structure 
of the and ^He gronnd states^^). Moreover, comparing the third and fonrth 
rows of the table, one can see a “nonlinear” effect of self-strengthening for the 3BF 
contribntion. In fact, the comparison of the resnlts presented in these rows of the 
table (see the second and the fifth columns) shows clearly that the binding energy 
is almost proportional to the weight Pin of the 6qN component in the total 3N 
wavefnnction. Thns, when the weight of the 6qN component increases, the 3BF 
contribntion, which is related directly only to this component of the total wave¬ 
fnnction, increases accordingly. However, the enhancement of the pnre attractive 
3BF contribution squeezes the 3N system and thus reduces its r.m.s. radius, i.e., 
the mean distance between nucleons, which, in turn, again increases the weight of 
the 6qN component. In other words, a some chain process which strengthens the 
attraction in the system arises. This process is balanced both by the weakening 
of the effective pairwise interaction dne to the dependence and by the repnisive 

It should be noted here that the relative contribution of the pairwise effective force W{E) to the 3N 

binding energy decreases noticeably when including 3BF (due to strengthening of the dependence 

arising from the pairwise forces). 
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effect of the orthogonalizing pseudopotentials included in each pair interaction. 

There are two another important stabilizing factors weakening the strong three- 
body attraction in the 3N system. First, the generation of the short-range repulsive 
vector cu-field, where all three nucleons are close to each other [25]. Since the 
cn-meson is heavy, this field is located in the deep overlap region of all three nucle¬ 
ons. In the present study, we omitted the three-body contribution of this repulsive 
cu-field. This repulsive contribution will keep the whole system from the further 
collapse due to the strong attractive 3N force induced by the scalar field. 

The second factor weakening slightly the effective 3N attraction is associated 
with the conservation of the number of scalar mesons generated in the 2N and 
3N interaction process. The problem is that TSE giving the 3BF contribution (see 
Fig. 4) arises due to the break of the a-meson loop, which induces the main 2N force. 
In other words, the a meson generated in the transition of pair nucleons from the 
NN phase state to the 6q state is absorbed either in the 6q bag with closing the loop 
or by the third nucleon, resulting in the 3BF contribution. Thus, the appearance of 
such a 3BF should weaken the attraction between nucleons in the pair. We carefully 
estimated the effect of the meson-number conservation for the TSE contribution on 
the total 3N binding energy. Its magnitude occurred to be rather moderate in the 
absolute energy scale (ca. 0.3-0.4 MeV), but quite noticeable within the whole TSE 
contribution. However, when the total nucleon density increases (and the relative 
TSE contribution also increases), the effect is enhanced. 

(v) Dependence of the two-body coupling constants A(e) upon the average momentum 
of other nucleon in 3N system (see, e.g., Eq. (29)) can be interpreted generally as 
a density dependence of the resulted many-body force in many-nucleon system. It 
is easy to show the appearance of the energy-dependent pairwise potentials of the 
above-mentioned type leads inavoidably to a repulsive many-body force. In other 
words, the effects of the two-body interactions of this type can be reinterpeted 
in terms of the conventional static interaction as additional contribution of the 
effective repulsive density-dependent many-body force. For example, if to remove 
the g^-dependence from the coupling constant X{E — q^/2m) of our two-body force 
(this g^-dependence leads to a weakening of the two-body force in a many-nucleon 
system, when g^ is rising), then the neglected g^-dependence must be compensated 
by an additional repulsive density-dependent effective three-body force. Thus we 
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can replace this energy-dependent two-body interaction by an effective static two- 
body potential (as it is usually done) plus a repulsive density-dependent 3BF. 

On the other hand, it is well known from the Skirme-model calculations of 
nuclei that just similar repulsive phenomenological density-dependent 3BF should 
be added to conventional 2N- and 3iV-forces to guarantee the saturation properties 
of heavy nuclei. Thus, in this respect the present force model is also in a qualitative 
agreement with phenomenological picture of nuclear interactions. 

9. CONCLUSION 

In this paper, we have developed a general formalism for the multicomponent descrip¬ 
tion of the three-body system with particles having inner degrees of freedom. We have 
applied our new approach to studying 3iV system with 2N and 3iV interactions based on 
the dressed dibaryon intermediate state and a-field generation. It has been shown that 
the DBM applied to the 3iV system results inevitably in new three-body scalar and also 
new (three-body) Coulomb forces due to the (strong -|- Coulomb) interaction between the 
dressed dibaryon and the third nucleon. These forces play a crucial role in the structure of 
few-nucleon systems and very likely in whole nuclear dynamics. Our accurate variational 
calculations have demonstrated that new 3BF gives a half of the 3iV binding energy, 
whereas the 3BF contribution in the traditional iVW-force approaches gives about 15% of 
the total binding energy. Thus, the suggested approach to the NN and 3iV interactions 
can lead to significant revision of relative contributions of two- and many-body forces in 
all nuclear systems. 

The developed model gives the precise value for the Coulomb displacement energy 
AUc of the A = 3 system. Two basic sources of this contribution, which differ from 
conventional force models, should be indicated here: 

the three-body Coulomb force between the dressed bag and the charged third nu¬ 
cleon; and 

quite significant correction to the kinetic energy of the system due to the np mass 
difference and high average kinetic energy. 

It should be emphasized that, contrary to other studies based on conventional force 
models (using the 2N and 3N forces generated via the meson-exchange mechanism), 
this explanation does not require any noticeable CSB effect, although our model is still 
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compatible with such effects. However, these CSB effects do not contribute remarkably 
to Ai?c in our approach. 

It is crucially important that the DBM leads to a significant non-nncleonic components 
in the 3N wavefnnction (8-11%), while this component in the denteron is ca. 3%, which 
resnlts in a reformnlation of many basic effects in few-nncleon systems and other nnclei 
as well. It is probable that the weight of such non-nucleonic components in heavy nuclei 
can be even higher with an increase in the mass number and nuclear density. 

There is a very specific new interplay between two- and three-body forces: the stronger 
the three-body force, the smaller the attractive contribntion of the two-body force to the 
nnclear binding energy! This gives a very important stabilization in nnclei and nnclear 
matter. By this way, a very natnral density dependence of nnclear interactions appears 
from the beginning. Thus, the general properties of the 3N system, where forces so much 
differ from any conventional force model, would appear also to be much differ from the 
predictions of any conventional model and, hence, from experiment. 

Therefore, it was very snrprising for ns to find that the static characteristics of the 
3N system in onr case tnrned ont to be very close to the predictions of the modern 
force model (snch as AV18 -I- UIX) and thns to experiment. This gives ns a good test 
of the self-consistency and accuracy of the new force model. However, predictions of 
the present NN- and 3iV-force model in other aspects will strongly be deviated from 
those for conventional models. First, these are the properties determined by the high- 
momentnm components of nnclear wavefnnctions. The point is that the system described 
by onr mnlticomponent wavefnnctions inclnding the dibaryon components explicitly can 
easily absorb high-momentnm transfers, which can hardly be absorbed by the system 
consisting of nucleons only. Therefore, to fit the experimental data corresponding to large- 
momentum transfers (~ 1 GeV/c), many types of meson-exchange and isobar currents are 
often introdnced into theoretical frameworks. However, these cnrrents are often nnrelated 
to the nnderlying force model. Hence, it is rather difficnlt to check the self-consistency of 
snch calcnlations, e.g., the validity of gange invariance, etc. 

Numerous modern experiments could corroborate these results. In particular, accord¬ 
ing to the recent experiments ^He(e, e'pp) [5] and their theoretical interpretation on the 
basis of fully realistic 3N calculations, the cross sections for the ^He(e, e'pp) process are 
nnderestimated by abont five times with a fnlly realistic 3N model and incorporation of 
final state interaction and meson-exchange cnrrents. This important conclnsion has been 
fnrther confirmed in recent experiments at the Jefferson Laboratory when the incident 
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electron beam energy has been increased up to E^. = 2.2 and 4.4 GeV [7]. The data of the 
two different experiments give a clear evidence of very strong short-range NN correlation 
in the ^He ground state. This correlation still cannot be explained within the traditional 
pattern for the 3iV system. 

In addition, our approach has recently been partially supported [67] from the other 
side by considering a model for 27r production in pp collisions at Ep = 750 and 900 MeV. 
The authors have found that almost all particle energy and angular correlations (e.g., 
TT+TT”, pp, npp, etc.) can be explained quantitatively by assuming that tt+tt” production 
occurs through the generation of an intermediate light a meson with the mass — 
380 MeV. These values generally agree with the parameters adopted in our NN model [8, 
24] and drastically disagree with the values assumed in OBE and other potential models. 

Very interesting general implication of the results presented here is their evident in¬ 
terrelation to the famous Walecka hadrodynamic model for nuclei [68]. It is well known 
that the Walecka model describes nuclei and nuclear matter in terms of the scalar a- and 
vector cu-fields, where the a-field gives the attractive contribution, while the vector cu-field 
balances this attraction by short-range repulsion. It is very important that both basic 
fields appear (in the model) as the explicit degrees of freedom (together with relativistic 
nucleons), in contrast to conventional meson-exchange models for nuclear forces, where 
mesons appear as the carriers of forces rather than as the explicit field degrees of free¬ 
dom. Our approach does include the a-meson (and potentially the cu-meson) degrees of 
freedom in an explicit form, similarly to the Walecka model. Moreover, since the average 
kinetic energy of the 3N system in our model is high (it is much higher than that in the 
conventional OBE approach), nucleon motion is closer to the relativistic case, and thus 
the similarity with the Walecka model gets even closer. 

There is also an additional strong argument in favor of a tight interrelation between 
our and the above Walecka-type nuclear model. Very recently, we have formulated [38] 
the dibaryon model for NN interaction in terms of relativistic effective field theory with 
the intermediate dibaryon being represented as a color quantum string with color quark 
clusters at its ends. This theory includes tt, a, p and cu-mesons as a dressing of the 
dibarion together with the iVA and AA loops. Thus, the 3iV scalar force introduced in 
the present work “by hands” can be derived in the field-theory Lagrange language within 
the effective field theory. Moreover, in the mean field approximation this effective field 
theory approach, being applied to nuclei, should result in the Walecka-Serot relativistic 
model with the dominating collective a-field, which couples the nucleons in a nucleus 
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together. 

Thus, the alternative description given here by the new force model looks to be more 
self-consistent and straightforward than the conventional OBE-type models. One aspect 
of this new pictnre is evident - the present model being applied to any electromagnetic 
process on nnclei leads antomatically to a consistent pictnre of the process as whole: 
single-nucleon currents at low-momentum transfers, meson-exchange currents (including 
new meson currents) at intermediate-momentum transfers, and quark counting rules at 
very-high-momentum transfers, because the model wavefunction includes explicitly mult¬ 
inucleon, meson-exchange, and mnltiqnark components. 

From all the above-mentioned argnments one can conclnde that the dibaryon concept 
of nnclear force advocated in the work resnlts in a new pictnre for nnclear strnctnre and 
dynamics. 
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APPENDIX 


OVERLAP FUNCTIONS BETWEEN 3iV SYMMETRIZED BASIS AND NN FORM 
FACTORS AND MATRIX ELEMENTS OF DBM INTERACTIONS FOR A 3iV 

SYSTEM 

1. The construction of basis functions 
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We use the following notations: rj(pj) is the relative coordinate (momentum) of the pair 
(jk), while Pj (qj) is the Jacobi coordinate (momentum) of the ith particle relative to the 
center of mass for the pair (jk), {i,j, k) = (1, 2, 3) or their cyclic permntations. Here the 
composite label 7 = {A / L Sjk S tjk\ represents the set of qnantnm nnmbers for the basis 
fnnctions: the angnlar momenta A and I correspond to the Jacobi coordinates r* and p^, 
respectively, Sjkitjk) is spin (isospin) of the two-body subsystem {jk), and L{S) is the 
total orbital momentum (spin) of the system. The normalizing coefficient in (A.3) is 






(A.4) 


7r(2A + 1)!!(2/ + 1)!! 

The spin-angular and isospin parts of the basis function are defined by a 
standard vector-conpling scheme: 


jrJMTT. = |{A,b ; T}{s,SA,(A,fc)s, : A} : JM), (A.5) 

= \tjtk{tjk)ti : TT^), (A.6) 

where Sj(= 1 / 2 ) and tj(= 1 / 2 ) are spin and isospin of the ith nncleon. 

Now we define the symmetrized basis functions as: 

= E (A-7) 

i=l,2,3 

SO that the total wavefnnction in an external (3iV) channel takes the form: 

= (A. 8 ) 

7,n 


2 . Nucleon-nucleon form factors 


The NN form factors in the separable DBM interaction and in the projectors 
9 ?^*^**'**'*"'(rj), corresponding to the orbital momentnm A*, spin Sd, the total angnlar mo¬ 
mentum and isospin (tdtdj of the subsystem (jk) (Ji = Ai -|- Sd), have the 

form: 


where 


^XiSd 


(d) = <^/(d) = 


DLri' exp(-i7^r2) Tj{ri) 


m 


(A. 9 ) 


f = {\i,Sd,JuMi,tddd,}, = \\iSd : JiMi), = \tjtk : tdtd,), (A.IO) 

and and rjm are linear and nonlinear parameters, respectively, of the Gaussian expan¬ 
sion. (In this Appendix we have altered the notation for the quantum numbers of the 
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NN form factor as compared with the main text of the paper: we have replaced Li —Aj 
for the orbital momentum and also included the isospin quantum numbers td,td^.) In 
the single-pole approximation the DBM includes only one form factor for each set /, so 
that index / determines the form factor uniquely. In the present version of DBM we use 
only OS', 2S, and 2D oscillator functions as the form factors. So, we need to expand in 
Gaussians the 2S function only. 


3. Overlap integrals 


The total overlap function (49) 


xLsD"' (“) ^ x/(») = (4’/(>)ii>r> 


(A.ll) 


772 


and also matrix elements (m.e.) of any DBM interaction include the overlap integrals 
between the form factors ipf and symmetrized basis functions 




This overlap integral consists of three terms: 


jfnn _ jfnn I r/,7« i jlin 

\i) (hi 


(A.12) 


(A.13) 


namely, one “diagonal” (/(pi) and two non-diagonal ones: 

m 

X [ Ti^ exp{-lrilrj)rjp\ exp(-a„r2 -(J^/(fi)|J7(rj,p.))d^D. 


Due to symmetry of the basis functions three overlap integras 


(A.14) 


U = 


1,2,3) are identical, so that further we present formulas for the case of i = 2. For 
example, 

x,(2) = ^CZ (ilpz + /|A,” + lim) (A.15) 
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(i) Diagonal overlap integrals I(2)2- 

^Y.<^22AeM-l3nPl) yf"’-'(P2)Ku, A, 

where 


^d^dz 


£ = (A.16) 


Jm 


0^2 - 


-*-1 A-I-.7/2 


2A+2 ^A+3/2 


I i J 
A S 31 Ji 
L S J 




(A.17) 
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[X] = 2X + 1, 

(A.18) 


1 2 

C^nm H“ 2^77*'’ 

(A.19) 



(A. 20 ) 


= (iA.i«2jrr.> |M 2 .>. 

(A.21) 

(ii) Non-diagonal overlap integrals: 


= ( 

_1)S4+S23 ^ G|iP*exp(-0i„rr7P2)3^/'^''^^(P2)^2l(fd,f23)A^2‘'*"A 

(A. 22 ) 

= 

(_1)A2+A Gl^plexp{-UJnmpl)yg'^^'^^{P2)^23{td,ti2)X2''^^% 

(A.23) 

where 

£ = {l,n,f,J,g,t,m}, 

(A.24) 


<^21 — ^*21(7) g, m, i -^^ 2 , L^,, /^4, j4, S'!, Qa}- 

(A.25) 

In Eq. (A.24) the summation is carrying out over all intermediate quantum numbers 
incorporated into the composite index Note that the overlap functions /( 2 )i and 1(2)3 
are distinguished by a phase factor and isospin functions only. 

The algebraic coefficients G 21 in (A.25) are equal to 


G 2 i( 7 . n. /, J, J, m, 0 = 4 “»):(Q„„) x 


J 9 .II 94 I VW|il|Sl|S23][SJ|J.]|A,-;]b4l[Jl,, 

^ r._Lr^_Lr._LQ ^ 


(2/i7 


Li -t-Z/3d-T4+3 


xr( 




((Ai-Li)Oj40|fi'o) 


(Aj — Li) Li Aj 

-S'o 3 gi 


X 


X 


23 


lie 
2 2^" 

i .S s. 


J (Aj —Li) 5(4 


J J L 


3a 


g (Aj-Li) gi g^ 



< 

1 

2 

> 


3a 

^4 


Li 

Sd 

9i 


L 

s 

J 

y 


In the formulas (A.22), (A.23), (A.26) the following notations are used: 


(A.26) 


Ainm — 




^nm — 


4/in 


(A.27) 


where 


13^ 3 1 73, 


(A.28) 
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The coefficients A in(A.26) are related to rotation of the basis functions from one 
Jacobi set to other one: 




Li L2 Ji \ I \ — Li I — L2 J2 


0 0 0 


X. 


'[A]![/]![A][/][L,][L2][A-L,][/-L2][ji]b2] 


[L,]\[L2]\[\-L,]\[1-L2]\ 

The rotation matrices P^m and Q„m in (A.26) have the forms: 


Li A-Li A 
L 2 I-L 2 I 
ji 32 L 


> . 


X 


0 0 0 


(A.29) 


‘Ifln 

0 1 


(A.30) 


^nm 


_ 1 
2 


2 

1 


2 2 

The overlaps between the basic isospin functions Tik are equal to: 






(A.31) 


= < 


\/( 2 «'t + l)( 2 i« + 1 ) 


1 1 r' 

2 2 jk 

h T tik 


X 


(-!)*«'= for (*j) = (13), (21), (32), (A.32) 

(_!)*'. for (^j) = (12),(23),(31). 


4. Conversion to momentum representation 


One of the advantages of Gaussian basis is the fact that Gaussian functions have the 
same form in both the coordinate and momentum representations. So expressions for the 
overlap integrals given below in the coordinate representation can be directly used for the 
calculation of m.e. of DBM interaction operators in the momentum representation. We 
use a “symmetrized” momentum representation: 

/(p) = j /(x)e--^. (A.33) 

Therefore, due to properties of the Gaussian functions, the (normalized) basis functions 
<h 7 n(pj, Qi) in momentum representation have the same form (A.3): 


(p„ q,) = NZpUI exp(-d„p2 _ ^^^2) WMTT. 


■ 7n \ 


Pi,qi 


ii) 


T 

tjk 


(AM) 
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where 


A(3n 

Moreover, as all the NN form factors (A.9) are the sums of Gaussians, the form of 
the overlap integrals (A.16)-(A.29) keeps invariable when passing from coordinate to 
momentum representation, if one replaces in these formulas: 


Qin ^ Oim /3n ^ Pni ^ • (A.36) 

Vm 

Below we use the symbols with tilde for designation of the corresponding quantities in 
momentum representation, e.g., = dn + etc. 


5. Matrix elements for DBM operators 


All quantities related to the non-nucleonic channels in DBM can be expressed in mo¬ 
mentum representation as the sums of integral operators with factorized kernels (see 
Eq. (36)): 

q^; E)ipf{pi). (A.37) 

Therefore the m.e. of such an operator is equal to the sum of the m.e. for one-particle 
operators qi; E) between the overlap functions X/(qj): 


M, = {4-„|Og5=“|>I>„> = ^ Mpp'iaj), (A.38) 

ff' 'yn,'y'n' ij=l,2^3 


fill 

where (*2j) are the corresponding basis m.e.: 




‘( 2 )i 


(A.39) 


Any scalar-isoscalar operator 0(q — q'), which does not depend on spin and isospin 
variables (e.g., the DBM two-body force, the projector, the 3BF due to a-exchange, the 
norm of non-nucleonic component), can be expanded into spherical harmonics as: 


0{q - q) = ^L{q, q)y^V(q')yLM(q) (A.40) 

LM 

In this case the spin-angular and isospin parts of the overlaps give: 




Sj'jSg'LSgL, 


^dz 


(A.41) 


(A.42) 
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Therefore, nine m.e.’s for such an operator M{i2j) = Mj 2 j {i,j = 1, 2, 3) can be reduced 
to radial integrals of four types (here we omit the indices /yn, f^'n! for brevity): 


M 122 = { — Ti2{t23,td)Rl22, 

M 322 = ( — 1 )^*^"^ T32{ti2Rd)Rl22, 

M 22 I = i — T2l{t23,td)R221, 

M 223 = { — T23{tl2-,td)R22l-, 

M 12 I = { — Tl2{t23-itd)T2l{t23,td)Rl21, 

M 323 = { — T32{t'i2Rd)T23{tl2,td)Rl21, 

Mi23 = i — ^)^‘^~^^^^^^"^^'Tl2it2^:td)T23{ti2,td)Rl21, 
M 32 I = (-l)^‘*+^^^+^^+^'r32(4,td)r2l(t23,td)i?121. 

Here, 


^121 - Y1 

G 2iG 21^J J'^gg'Rg i ^nrm 0^, 

(A.43) 

J’ J’' gg' ,mm' 

'tt' 


R 122 = 

G2lGl2^JJ'Sg/lRl\uJn'm',Pn, 0), 

(A.44) 

J J' g' ,mm't' 


R 22 I — 

G 22 G21^JJ' 5 gl'Rl'vi.Pn'^^nm^ 0), 

(A.45) 


J J'g,mm't 


R222 — ^ G22^22^JJ'^WPn', f^n'iO), 

J J' ,mm' 


(A.46) 



R^^{u\u-0)= / 


(A.47) 


0 0 


Below we give the explicit formulas for the radial integrals R^^ for all specific terms of 
the DBM interaction. 

Projector. The total projector onto the state ip'l", has the form: 




Ji^i^d^dz 

Sd 


(A.48) 


After expanding the ^-function into partial waves: 


<^(q-qO = '^yiMi^)yLM{(i)S{q - q')/q‘^, 

LM 


(A.49) 
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the corresponding operator Ol in Eq. (A.47) is reduced to 

and the radial integral for projector takes the form: 

r(^^) 




2(0.' + ^)*^’ 


(A.50) 


(A.51) 


where r(a;) is the gamma-function. 

Effective two-body DBM interaction. According to Eq. (29), the two-body DBM 
interaction in the 3N system (between nucleons 1 and 3) has the form: 


where 


therefore 


1E2= V lEi,(p',P2;q2,q2;i?), 

f * i * 


= E - q')AiA,(£-9V(2A))(43j";|, 


Mi 






(A.52) 


(A.53) 


(A.54) 


L 

In the present version of DBM we employed a rational approximation for the energy 
dependence of the coupling constant X'h^ilE) [24]: 

Eq + aE 


\iAE) = Aii,(0). 


Eq — E 


(A.55) 


where the parameters Eq and a can be taken to be the same for all A’s. We found that 
this simple rational form can reproduce quite accurately the exact energy dependence of 
the coupling constants A;^^,(E') calculated from the loop diagram in Eig. 1. Therefore, in 
the 3N system the corresponding coupling constants Ay^ take the form: 


“ -^A'Ai(O) -a + (a + 1) 


Eq 


Eq — E 1 -|- 


(A.56) 


2fh(Eo—E) , 

The first term in Eq. (A.56) leads to the expression for the radial matrix element like 
(A.51). Eor calculating the second term we expand the function 1/(1 + q"^) into a sum of 
Gaussians: 

(A.57) 


M 


1 + gVgg 

where gg = 2m{EQ — E) > 0 (for E < Eq) and the expansion parameters are 

universal constants. Then the total m.e. for two-body DBM interaction kE 2 takes the 
form: 
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+ (a + 1) 


En 




r( 






(A.58) 


Norm of 6 qN component. The norm of 6 qN component of the total 3N wavefnnc- 
tion determined by Eq. (34) can be expressed via a sum of the m.e. of the operator: 

A/'(q', q) = (5(q' - q). (A.59) 

For the energy dependence such as in Eq. (A.55) the derivative takes the form: 


Therefore, 


r( 




'M 


M'M 


2{uj'+ uj + {6 m + (^M)/qQ) 2 

(A.61) 


Three-body force due to OME. 

When calculating the matrix elements for 3BF due to OME (Eq. (37)), viz.: 


OMEij/J'iJirM 


lf^;:;(q',qdi?)= / dk- 


Biik') 




E - E. 

2m 


Bt:(k) 






(A.62) 


we use a similar trick as in the calculation of the norm for the QqN component. It enables 
us to exclude the vertex functions Bj^,{k) from the formulas for the matrix elements. By 
replacing the product of propagators in the integral (over the meson momentum k in 
Eq. (A.62)) with their difference, one obtains the following expression free of the vertex 
functions: 




(E-e(k)-f-)(E-e(V)-^) 


dk = 




g/2 _ ^2 


= A\i:,,^{q',q) (A.63) 

This quantity is the finite-difference analogue of the derivative of A with respect to g^, 
which, in the case of the energy dependence (A.55), takes the form: 

1 1 


A\{q , q) — A(0) Eo(l-l-a)- 


(A.64) 


E — q^/2m E—q'"^/2m 

Thus, the matrix elements for OME can be found without explicit usage of the vertex 
functions Bj/,{k). 
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Three-body force due to OSE. The exchange operator for scalar mesons does not 
include any spin-isospin variables. Therefore, Eq.(A.62) can be simplified and, in view of 
the energy dependence given in Eq. (A.55), reduces to the form 


E) = £„(! + a)x 


1 


X 




Using the expansion of the OME interaction into partial waves: 

T ,— W ^ 

(q) - + ml 2qq' 

where is mass of a meson, Ql{z) is Legendre function of the second kind, and 


(A.65) 


(A.66) 


2 : = 


+ q^ + ml 

2q'q 


one gets the following radial integral for OSE m.e.: 


OSEttx/'/n_2 \Ji + 




(E - Eo)2 


'y ) l3j^/l3jiiTZ(t ,uj -\- 9 ,t, Co -\- 9] L, 


m„ 


M'M 


(A.67) 


where 


7l{t',u L,m) = f f{qY^^e '^‘^dqdq (A.68) 

J J 

We calculate integrals like those in Eq. (A.68) in the following way. In the integral 
(A.68) t' > L, t > L and it can be shown that t' — L and t — L are the even numbers. 
Introducing the auxiliary indices k and k\ so that 

t = L + 2k t' = L + 2U, 

the integral (A.68) can be written as 

( Q \k' / Q \ ^ 

— (~^) E.(L, cu', L, cu; L, m). (A.69) 

The last integral with t = t' = L is easily calculated in coordinate representation. Using 
the well-known formulas 


Ql{ 


q'^+q’^+m? \ °° 


^—mp 


-= / jdq'p) - - jL{qp)p^dp 

2qq' J p 


(A.70) 


and 


00 

1 j = 


(A.71) 


0 
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we get the result: 


7^(L, o;', L, o;; L, m) = — 

where 

OO 

_ j p2L+2^^^-bp^^p_ 
0 

This integral of Yukawa potential is calculated by recursions: 

1 771 


where 


Zq{x) = — erf(a;)), 

Zi{x) = 2x^(1 - Zo{x)), 

Zl{x) = 2x^ {{k - 2 )Zl-2{x) - Zl-i{x)). 


(A.72) 


(A.73) 


(A.74) 


(A.75) 

(A.76) 

(A.77) 


These recursions (especially for large L) for functions Z{x) and also expressions (A.75)- 
(A.76) for Zq and Zi cannot be used for large values of x. So, at large x we use the 
following asymptotic series: 




2by {2L + 2i - 1)\\{L + i)\ 


(ATS) 




26 y (2L + 2i + l)!!(L + i)! 


(A.79) 


Three-body force due to OPE. For OPE, we take the interaction operator in the 
standard form 




(A.80) 


where and t^') are the spin and isospin operators of the third (Ah) nucleon, whereas 
Sd and Trf are the operators of the total spin and isospin of the 6g bag, respectively. We 
found that the contribution of OPE is so small that it is sufficient to include only S waves 
in its evaluation. In this case, the central part of the OPE interaction only is remained: 

(A.81) 

The spin-isospin matrix element is nonzero only for a singlet-triplet transition: 

(Ad = 0,Trf = l|i(a« ■ Sd)(T« ■Td)\Sd = l,Td = 0) = ^. (A.82) 

o y 
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Then, the matrix element of the OPE contribution for 5'-waves takes the form (for Sd = 
0,S'^ = 1 or S'd = l, S'(; = 0): 

4 




X 


^ ^ 3j^iB_x4TZ{0^ uj + 0,0, o) + 0, 




(A.83) 


M'M 


Here, we take the vertex functions Bq and Bq differing from each other by a constant 
only. Therefore, using Eq. (A.63), one can exclude these functions from the formula for 
the matrix element. 

Three-body Coulomb force. The m.e. of the operator for three-body Coulomb 
force (60) (for point-like charges) can be expressed in terms of the integrals over the 
overlap functions X/(q): 


' T (B-£„-£)(q-qOM£-J5o-i2 

(A.84) 

where (l + Tg j/(2) is the operator of the nucleon charge (nucleon with number 2). The 
isospin part of this m.e. is equal to: 


( 2 ) 


72 - dqdq^, 


T 


Coul 


iu) = 

tdz 


l-kTg 


( 2 ) 


(1 + tdj 


A*'* > = 


1, for td = 0, 


(A.85) 


i, for td = 1. 

Thus, for calculation of 3BE Coulomb m.e. one can apply the formulas Eqs. (A.69)-(A.79) 
for isoscalar operator with this additional isospin factor (A.85): 


CouItt//'/\ _ r , ,^„2 \ Ji /n\ -^o(l + «) 


Co; ^) = 5,2,, (0) 


{E — Eq 


B'j^,BM'J^'^°"\t\Co'+9\t,Co+9; L). 

(A.86) 


M'M 


Here the Coulomb integral ,t,u; L) for the point-like charges is obtained from 

Yukawa integral IZ (A.68) by putting m = 0: 


77^°"^(T, uo\ t, uo; L) = n{t', u', t, u; L, 0). (A.87) 

Hence these Coulomb integrals are reduced by differentiating (see Eq. (A.69)) to the 
integrals: 


n^°'^\L,u\L,u-,L) = 


- f p2L+2ig-4(a; 


{Auouo'Y+^B 


jC-^((o;-i+o;'-^)/4), 


(A.88) 
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JL°'‘'(b) = (A-89) 

Now, we can replace the Coulomb potential 1/p between the point-like charges in the 
integrand in Eq. (A.88) with the corresponding Coulomb potential between the “smeared” 
charges: 


^ y p2L+2 
0 

The latter integral is evaluated analytically in the form of a finite sum: 


(A.90) 




1 1 




k\ 


{2L-2k + l)\\ 


2 fli+l ^ ^ (^k'jk+l 2L-k(^b _|_ 


(A.91) 


where are the binomial coefficients. 
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